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Abstract 


Analytical  and  computational  models  of  fiber  reinforced  composite 
materials  are  constructed  using  a  mathematical  procedure  called  ho¬ 
mogenization.  The  procedure  is  systematic;  i.e.,  one  can  estimate  the 
precision  of  the  approximations  to  the  field  equations  of  the  compos¬ 
ite.  Moreover,  it  retains  interaction  effects  due  to  the  microstructure 
of  the  composite  in  the  macroscopic  approximations.  A  software  sys¬ 
tem,  MeMCAP,  based  on  finite  element  methods  implements  the  method- 
ok  ^y  for  evaluation  of  macroscopic  eSectiv.^  moduli  of  the  composite 
and  for  computation  of  the  microscopic  stress  and  strain  fields  acting 
at  the  fiber-matrix  interfaces.  The  software  system  is  designed  for  easy 
use  by  engineers  who  need  not  be  familiar  with  the  underlying  analyt¬ 
ical  techniques.  Interaction  takes  place  through  a  menu-driven  control 
structure  which  requires  specification  of  the  material  properties  of  fiber 
and  matrix  and  of  the  “geometry”  of  a  typical  cell  of  the  composite. 
The  homogenization  method  produces  an  effective  parameter  model  of 
the  macroscopic  behavior  of  the  material  (longitudinal  and  transverse 
Young’s  moduli,  Poisson  coefficients,  etc.),  and  a  description  of  the  mi¬ 
croscopic  distributions  of  stress  within  the  cell  -  especially  at  the  fiber- 
matrix  interface.  The  system  can  treat  composites  with  various  fiber 
'  shapes  and  packing  arrangements.  It  can  also  treat  multi-ply  laminates. 
The  program,  written  in  Fortran,  is  fast  and  efficient.  Extensions  of 
the  work  to  encompass  dynamical  phenomena  such  as  the  propagation 
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roperties,  and  estimation  of  the 
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Introduction  and  Project  Summary 


The  objective  of  this  project  is  the  development  of  analytical  methodology  im¬ 
plemented  in  a  software  tool  for  the  modeling  and  design  of  composite  materiads. 
Otit  ultimate  objective  is  to  provide  a  software  system  for  engineering  level  eval¬ 
uation  of  composites  in  a  form  which  is  easy  to  use,  reasonably  comprehensive, 
flexible,  extensible,  portable  and  computationally  efficient.  To  this  end  we  have 
developed  a  prototype  program  for  the  evaluation  of  effective  elasticity  moduli  for 
flber-reinforced  composites  which  includes  the  capability  to  compute  an  approxi¬ 
mation  to  the  microscopic  distribution  of  stresses  at  the  flber-matrix  interface. 


1.1  Modeling  Composite  Materials 

Direct  computation  of  the  detailed  behavior  of  composites  imder  loads  is  virtu¬ 
ally  impossible  using  conventional  methods,  including  sophisticated  finite  element 
codes.  This  is  a  consequence  of  the  large  number  of  degrees  of  freedom  in  the 
heterogeneous  material.  A  computational  procedure  which  has  proven  useful  is  to 
derive  a  model  representation  with  material  coefficients  which  do  not  vary  (rapidly) 
in  spatial  variables  -  in  effect  to  find  a  homogeneous  material  whose  macroscopic 
behavior  approximates  that  of  the  composite  in  a  specifc  context.  Many  techniques 
have  been  used  for  the  derivation  of  homogeneous,  continuum  approximations  of 
composite  materials.^  In  this  report  we  shall  focus  on  the  homogenization  method, 
particularly  as  developed  in  [5,6,7]  and  in  the  general  sources  [12,35,65].  This  is 
a  powerful  mathematical  technique  for  the  analysis  of  physical  systems  in  which 
there  are  two  or  more  scales  upon  which  spatial  or  temporal  (or  both)  phenomena 
occur.  In  a  composite  medium  the  natural  scales  are  the  small  (microscopic)  scale 
of  interfiber  spacing  in  fiber  reinforced  materials,  or  the  mean  free  path  length  be¬ 
tween  particles  in  particle  reinforced  compc^ites;  and  the  large  (macroscopic)  scale 
characterizing  the  overall  dimensions  of  the  structure  formed  by  the  medium. 

This  report  presents  two  basic  results  from  application  of  the  homogenization 
method  to  the  analysis  of  fiber — ^reinforced  composite  materials: 


^See  among  many  other  references  (1,2,3,22,23,74). 
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1.  When  the  period  (dimension)  of  a  basic  “cell”  of  the  periodic  structure  ap¬ 
proaches  zero,  the  fields  of  deformations  and  stress  of  the  composite  tend  to 
those  corresponding  to  a  homogeneous  (anisotropic)  structure.  The  complete 
set  of  moduli  for  the  homogenized  structure  can  be  computed  in  terms  of  the 
elastic  moduli  of  the  constituents  (fibers  and  matrix)  and  the  parameters  de¬ 
scribing  the  geometrical  layout  of  a  single  “period”  or  “cell”  of  the  composite 
structure. 

2.  By  retaining  additional  terms  in  the  homogenization  asymptotic  analysis,  we 
can  develop  a  picture  of  the  local  distribution  of  stresses  in  the  material, 
e.g.,  at  the  interfaces  between  the  fibers  and  matrix.  (See  Figure  1  and  the 
examples  in  Section  4.) 

The  homogenization  procedure  is  syatem<Uie;  i.e.,  one  can  estimate  the  precision 
of  the  approximations  to  the  field  equations  of  the  composite.  Moreover,  unlike  the 
“rule  of  mixtures,”  it  retains  interaction  effects  due  to  the  microstructure  of  the 
composite  in  the  macroscopic  approxhnations.  Numerical  evaluation  of  both  the 
homogenized  moduli  and  the  distributions  of  local  stresses  csm  be  based  on  finite 
element  analysb  of  elliptic  partial  differential  equations. 

The  homogenization  procedure  is  general  in  the  sense  that  it  can  treat  both 
materials  with  a  regular,  periodic  infrastructure  and  materials  with  a  random  in- 
frastructxire.  The  general  form  of  the  analysis  is  similar  in  both  cases,  as  is  the  form 
of  the  expressions  for  the  effective  parameters  in  the  approximations;  of  course,  the 
details  of  the  derivations  are  quite  different.* 

In  this  Phase  I  project  we  have  developed  a  prototype  software  system  based 
on  finite  element  methods  to  implement  the  homogenization  methodology  the  eval¬ 
uation  of  macroscopic  effective  moduli  of  composites  with  periodic  infrastructure 
aind  for  computation  of  the  microscopic  stress  and  strain  fields  acting  at  the  fiber- 
matrix  interfaces.  The  software  system  is  designed  for  easy  use  by  engineers  who 
need  not  be  familiar  with  the  \mder lying  analytical  techniques.  Interaction  takes 
place  through  a  menu— driven  “supervisor”  program  which  requires  specification  of 
the  material  properties  of  fiber  and  matrix  and  of  the  “geometry”  of  a  typical  cell  of 

*Tlie  general  theories  for  the  periodic  and  random  cases  are  developed  in  [12]  and  |69j,  respec¬ 
tively. 
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Figure  1:  Microscopic  stress  field  in  a  cross  section  of  a  typical  cell  computed  by 
homogenization. 
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the  composite.  The  homogenization  method  produces  an  effective  parameter  model 
of  the  macroscopic  behavior  of  the  material  (longitudinal  and  transverse  Young’s 
moduli,  Poisson  coefficients,  etc.),  and  a  description  of  the  microscopic  distributions 
of  stress  within  the  cell  -  especially  at  the  fiber-matrix  interface.  The  system  can 
treat  composites  with  various  fiber  shapes  and  packing  arrangements.  It  can  abo 
treat  multi-ply  laminates.  The  program,  written  in  Fortran,  b  fast,  efficient,  and 
portable.  It  can  be  readily  extended  to  treat  other  types  of  composites,  and  it  can 
be  enhanced  to  incorporate  a  database  of  composite  properties  for  comparbon  with 
predictions  of  the  analytical  modeb.  We  plan  to  implement  these  extensions  and 
enhancements  in  the  second  phase  of  thb  project. 


1.2  Project  Summary 

The  results  of  the  project  fall  into  two  classes: 


1.  Analytical  results  on  the  derivation  of  the  "homogenized  approximations” 
for  the  elastic  response  of  fiber-reinforced  composites  in  the  periodic  case, 
including  asymptotic  expansions  of  the  stress,  strain,  and  dbplacement  fields; 
and 

2.  Implementation  of  the  effective  parameter  modeb  and  field  expansions  in  a 
menu-driven  software  system. 


The  analytical  results  are  summarbed  in  the  next  section.  The  softwue  system 
b  described  in  section  3.  Sample  numerical  results  from  the  software  system  are 
given  in  section  4  and  Lutings  1  and  2  at  the  end  of  the  report.  Recommendations 
for  further  research  are  given  section  5. 


1.3  Comments  on  Some  Related  Work 


The  need  for  toob  for  stress  analysb  of  laminates  formed  from  composites  as  a 
basb  for  the  evaluation  of  strength  and  failure  properties  has  been  emphasized  by 


I 
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many  researchers,  see  for  example  [71].  As  pointed  out  in  [71],  difficulties  have  been 
encountered  in  treating  systems  containing  multiple  layers.  As  we  show  in  section 
3  (see  also  [5]),  the  homogenization  method  readily  handles  such  cases,  even  when 
the  number  of  layers  in  the  laminate  is  large. 


2  The  Homogenization  Method  for  Modeling  Com¬ 
posite  Materials 

2.1  Formulation  of  the  Model 


Consider  the  problem  of  characterizing  small  scale  elastic  deformations  of  a  com¬ 
posite  material  which,  in  its  undeformed  state,  occupies  a  region  H,  open  and 
bounded  in  with  regular  boundary  T.  Suppose  the  boundary  b  partitioned 
into  r  =  r©  U  Ti.  (See  Figure  2.)  Suppose  the  material  is  constrained  (not 
to  move)  along  the  F©  portion  of  the  boundary,  and  that  a  surface  force  g  = 
[gi(zi,X3,X3),g2(2ii23t2:3),g3(xi,X3,Z3)]  b  incident  on  the  material  along  the  Fi 
portion  of  the  boundary.  Suppose  in  addition  there  are  volumetric  forces  /  = 
[/i(®i»*2»*3)>/2(®i>*2i*3)i/3(®ija:jia:3)]  acting  on  the  interior  of  the  body. 


Let  u(xi,X2,X3)  =  [u2(xi,X2,X3),tt2(xi,X2,X3),«3(a:i,X2,X3)]  be  the  dbplace- 
ment  field  of  the  body  subjected  to  the  force  densities  f,g.  Then 


1  <  »,J  <  3 


(1) 


defines  the  strain  tensor.  The  stress  tensor  b  a,-y  which  satbfies  the  equilibrium 
equations 


^^'+/.=Omn 

(2) 

Oijnj  =  gi  on  Fi 

(3) 
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Figure  2:  Elastic  body  subjected  to  forces. 

Stress  and  strain  axe  related  by  the  pointwise  constitutive  relation® 

<T,-,(u)  =  a,-,«(xi,i2,23)eM(u),  1  <  *,y  <  3  (4) 

The  components  of  the  elasticity  tensor  aiju{xi,X2,X3)  axe  bounded  and  satisfy 
symmetry  a,,fc«  =  a^iu  =  Oijtk ,  l<iyJ,k,t<Z 

positivity  >  a  X) 

for  some  a  >  0,  and  for  all  such  that  =  £y,- 

In  composite  materials  with  a  periodic  infrastructure  of  fibers,  cells,  or  layers, 
etc.,  the  elasticity  tensor  will  be  a  periodic  function  of  the  spatial  variables.  The 
homogenization  method  is  suitable  for  treating  the  case  when  the  infrastructure  of 
the  composite  is  “fine”  relative  to  the  macroscopic  dimensions  of  the  structure.  The 
periodic  structures  treated  here  can  be  simple  identical  rectangular  or  hexagonal 
arrangements,  or  they  can  be  more  complex  shapes  as  shown  in  Figure  3.  The 
only  constraint  is  that  the  opposing  faces  of  a  typical  “cell”  of  the  structure  which 
correspond  through  a  translation  can  be  identified  in  pairs. 

*In  the  following,  we  nee  the  convention  that  repeated  indices  are  summed. 
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Figure  3;  Typical  periodic  cross-sections  of  composites  for  the  homogenization 
method. 

Let  y  be  a  dimension  of  a  typical  cell  of  the  structure,  and  suppose  that  the 
structure,  i.e.,  the  overall  composite,  consists  of  a  large  number  of  cells.  Let  c  be 
a  small  dimensionless  parameter  that  b  the  (homothetic)  ratio  of  y  to  a  typical 
dimension  of  the  structure  as  a  whole.'*  If  the  spatial  variable  x  €  b  used  to 
describe  the  macroscopic  variations  in  the  material,  then  we  use  y  =  x/e  as  the 
microscopic  spatial  scale  which  describes  rapid  variations  of  the  elasticity  across 
the  structure.®  Hence,  it  b  appropriate  to  identify 

X 

0’ijkl{^U^2,X3)  =  OiMv)^  ^  7' 

Introducing  the  shorthand  notation, 

X 

«(S/)  =  tt*(a:i,X2,X3)  =  a{-) 

*For  example,  if  the  composite  structure  is  a  fiber  reinforced  rod  of  length  L  and  thickness  W, 
and  if  a  typical  cell  of  the  rod  is  defined  to  be  a  section  of  the  rod  including  a  single  fiber  and  a 
portion  of  the  surrounding  matrix,  and  if  the  cell  cross  section  as  length  Y,  then  c  =  Y/W.  Hence, 
c  is  dimensionless,  and  assuming  there  are  many  cells  in  the  cross  section  of  the  rod,  0  <  c  <<  1. 

•Using  the  composite  rod  example,  if  we  move  the  small  distance  of  x  =  «  units  through  the 
cross  section  of  the  rod,  we  will  pass  from  matrix  to  fiber  and  back  to  matrix,  incurring  a  significant 
variation  of  elasticity.  This  corresponds  to  a  movement  of  one  unit  in  the  fast  spatial  scale  y. 
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a  =  {dj},  e(tt)  =  {ew(tt)} 
the  constitutive  law  (4)  is  written 

=  o-iiki{^i,X2yX3)eu  (6) 

where 

<(«)  =  K  W).  «</(«)  =  I  (^  +  ^) 

When  an  ambiguity  is  possible,  either  ex(u)  or  ey{u)  will  be  specified  depending  on 
whether  the  derivative  is  with  respect  to  x  or  y.  The  boundary  conditions  are: 

tt  =  0  on  To  djnj  =  y,-  on  Fi  (7) 

The  problem  posed  by  (l)-(4)  has  a  tmique  solution  which  depends  on  e  and 
which  we  shall  designate  as  u*;  a  stress  field  a*  may  be  associated  with  this: 

<r‘  =  tt‘(x)e(tt‘)  (8) 

When  e  is  small,  it  is  very  difficult  to  calculate  numerically  since  there  are  a 
large  number  of  heterogeneities  in  the  elastic  medium.  Therefore,  one  usually  tries 
to  obtain  an  expansion  of  the  solution  u*  and  the  stress  field  o*. 

2.2  The  homogenization  method 

The  homogenization  method  results  in  replacing  the  pointwise  constitutive  rela¬ 
tion  (6)  by  a  linear  relation  (with  coefficients  constant  in  space)  between  the  "mean 
values”  of  the  stress  and  strain  tensors.* 

*It  is  important  to  emphasise  that  the  method  applies  not  only  to  systems  with  a  periodic  struc¬ 
ture,  but  also  to  systems  with  a  structure  which  varies  randomly  in  space.  The  random  variations 
must  be  stationary  in  a  strong  sense,  and  the  detuls  of  the  arguments  differ  from  the  periodic  case; 
however,  the  principles  are  similar.  See  [59]  for  the  general  theory  in  the  random  case  and  |10,50] 
for  extensive  applications. 
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2.2.1  Asymptotic  Expansions 
The  solution  depends  on  two  scales: 


(i)  The  first  b  the  macroscopic  dimension  (scale)  of  fl  which  affects  the  solution 

through  the  forces  applied  to  the  body  as  whole  and  the  conditions  at  the 
boundaries. 

(ii)  The  second  b  due  to  the  period  of  the  microscopic  system  infrastructure  which 
characterbes  the  local  internal  forces  and  stress  dbtributions. 


Thb  jiistifies  looking  for  an  asymptotic  expansion  of  the  solution  the  form: 


u*  =  y)  +  6tt'(i,y)  +  c’tt’(i,  y)  +  •  •  •  (9) 

where  the  u*(x,y)  are,  for  each  x  €  fl,  IT-periodic  functions  with  respect  to  the 
variable  y  €  Y.  Then  y  =  x/c  in  (9)  b  the  “fast”  local  scale  of  the  internal 
structural  variations,  i.e.,  x  scaled  by  the  (normalbed)  dbtance  between  fibers.  An 
expansion  of  the  deformation  field  e(u*)  may  also  be  identified 

c(u‘)  =  ie„(tt®)+e(u®)+e,(tt‘)  +  €[e,(u^)  +  ej(tt*)]  +  -..  (10) 

and  the  associated  stresses  field  has  the  representation 

O*  =  y)  +  y)  +  ca*  H -  (11) 

with 

a®  =  a(y)e„(tt®) 

c'  =  a(y)[ev(u‘)  +  e,(u®)] 

<^^{x,y)  =  [e»(u*)  +  e,(u*)] 


The  equilibrium  equations  (2)  applied  to  a*  give 
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or  in  a  more  condensed  form 

div  <T*  +  /  =  0 

Given  the  expansion  (11)  of  a*  we  have 

-^divyo“  +  -  (div„<T^  +  div*a®)  +  div^o*  +  div*<7'  +  /  H —  •  =  0 

=  =  {M} 

The  boimdary  conditions  (3)  are  treated  in  the  same  way: 

•  n  +  •  n  —  +  ccr*  •  n  +  . . .  =  0,  x  6  Ti,  y  EY 

Finally,  the  conditions  (7)  mean  that 

u®  +  €tt^  +  6*u*  + . . .  =  0,  X  €  To,  y  EY 

By  making  the  variotis  powers  of  e  zero  we  obtain: 


( 


divyO®  = 


o®  =  o(y)e„(u°) 


{ 


I 


divyO^  +  div,a°  = 


=  «(y)[ev(“^)  +  e.(tt®)] 


divy<7^  +  diVjO-^  +  /  = 


=  a(y)[ey(tt*)  +  e.(tti)] 


The  equations  (14)  and  (15)  will  be  used  later. 


(12) 

(13) 


(14) 

(16) 

(16) 

(17) 

(18) 


2.2.2  Analysis  of  the  Asymptotic  Expansion 

The  systems  (16)-(18)  contain  differential  operators  in  y.  They  are  equations  with 
partial  derivatives  relative  to  the  macroscopic  spatial  scale  x  with  the  associated 
period  Y;  and  so,  the  unknown  factors  are  y-periodic  functions. 
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System  (16):  This  leads  immediately  to: 

=  0,  ft®  =  tt®(x) 

System  (17):  In  view  of  (id)  it  reduces  to: 

divyff^  =0,  =  a(y)(ey(u^)  +  e*(u®)] 


(19) 


(20) 


The  deformation  es(u®)  is  a  function  of  x  and  not  y;  therefore,  it  plays  the  role 
of  a  parameter  with  respect  to  the  differential  system  in  y.  As  a  consequence  of 
linearity,  and  may  therefore  be  written  in  the  form: 


I 


=  «“(y)eM(u®) 

=  x“(y)eM(t‘®) 


(21) 


where 


divys***  =  0 


=  tt(y)(»’“  +  «»(x“)] 

Mu")  =  |(^  +  IS) 

is  y-periodic 


(22) 


The  tensor  r**  is  the  unit  tensor. 


It  can  be  shown  that  the  system  (22)  determines  the  vector  X^(y)  to  within  an 
additive  constant. 

For  any  function  ^  =  ^(x,y),  we  define  the  averaging  operation 

The  solution  of  (17)  is  given  by, 

<f'{x,y)  =  o(y)[r“  -  e»(x“))eM(tt“)  (23) 

and  taking  the  mean  value,  we  obtain, 

<  <’«  >=  Stf'uK) 


(24) 
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where 

?“  =<  «(«(»)  >  -  <  oi/mM'mCx")  > 


(25) 


System  (18):  It  suffices  to  take  the  mean  over  a  period  Y  in  the  first  equation 
to  obtain: 

div*  <  <7^  >  +/  =  0  in  n  (26) 


If  we  introduce  E  =<  >,  we  have 

(  divxE  +  /  =  0  in  n 

1  Si/  =  9^^eM(u®) 

Using  equation  (15)  and  taking  the  mean  on  V  in  (14),  we  obtain: 


(27) 


f  tt®  =  0  on  To 
(  E  •  n  =  y  on  F  f 


(28) 


The  system  (27)  with  boundary  conditions  (28)  is  a  well  posed  elasticity  problem; 
the  equilibrium  equations  are  unchanged,  as  well  as  the  boundary  conditions.  The 
elastic  constitutive  relation  is 


Si/  =  9^*e«(u®)  (29) 

It  is  homogeneoiis  since  the  coefficients  given  by  (25)  are  independent  of  z  G  H. 
These  coefficients  define  the  equivalent  "homogeneous”  material.  They  are  called 
homogenized  or  effective  eoeffieiente. 

Note:  The  effective  parameters  q  (25)  include  not  only  the  averages  of  the  a{y), 
the  “rule  of  mixtures"  which  one  expects,  but  also  terms  which  “correct"  for  local 
interactions.  These  correctors  are  usually  omitted  in  naive  effective  medium  theo¬ 
ries  fsj. 

The  stress  field  E  =  {I^t/}  is  called  the  maeroseopie  stress  field.  It  is  defined  by 
E  =<  >.  The  strain  field  E  =  Cs(u®)  is  called  the  maeroseopie  strain  field  and 

satisfies 

E  =<  e,(tt®)  +  e^(u^)  > 
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It  can  be  proved  that  the  homogenized  coefficients  qff  satisfy 
/  Qii  =  Qkei=  Qiikt) 

1  >  otiSijSij,  Vsy  =  sy,-  some  Oj  >  0 

This  shows  that  are  reasonable  elastic  coefficients  and  that  the  macroscopic 
scale  problem  (27)  (28)  has  a  imique  solution. 

This  homogenized  system  describes  the  macroscopic  deformation  and  stress 
fields  of  the  system  subject  to  external  forces  and  bormdary  conditions.  It  is  the 
effective  parameter  representation  of  the  medium.  Since  we  are  interested  in  the 
microscopic  stresses  which  characterize  the  fiber  matrix  interface,  we  must  push  the 
analysis  a  little  further  to  reveal  this  information. 


2.2.3  Microscopic  Fields  -  Localization 

The  stress  field  <r^(x,y)  is  the  first  term  of  the  asymptotic  expansion  (11)  of  the 
stress  field  c*{x)  solution  of  the  initial  exact  problem.  The  field  a^{x,y)  is  called 
the  microscopic  stress  field.  If  we  imagine  that  at  each  point  x  6  H,  there  is  a  small 
eY  period  with  its  composite  structure,  i.e.,  a  “cell”  of  the  overall  structure,  then 
gives,  for  x  kept  fixed  in  O,  a  stress  field  in  this  period. 

It  can  be  shown  that  c*{x)  —  <t^(x,x/6)  tends  to  zero  when  c  tends  to  zero.^ 
This  proves  that  <r‘(x,x/«)  is  a  good  approximation  to  ct*(x)  when  c  is  small.  The 
microscopic  stress  field  c^{x,y),y  =  x/«  can  be  calculated  as  follows: 

(i)  First,  we  obtain  the  six  x*^(y)  vector  fields  on  Y,  each  one  associated  with 

tensor  r**  =  r“.  These  six  vector-fields  are  solutions  of  problem  (22)  which 
is  an  elasticity  problem  on  the  large  scale  (period)  Y. 

(ii)  From  the  vector  fields  x*^  we  get  the  homogenized  coefficients  by  for¬ 

mula  (25). 


’In  the  (D)  norm;  that  is,  liin(_o  /„  \<r*{x)  —  a’(i,  x/c)|(Iz  =  0. 
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(iii)  We  solve  the  macroscopic  scale,  homogenized  elastic  problem  (27)  (28)  on  Q 

for  the  macroscopic  stress  field  £(x)  and  the  macroscopic  strain  field  = 

E{x)^  for  X  6  n. 

(iv)  The  localization  procedure:  Using  formula  (23),  we  can  calculate  a^(x,y). 
For  X  fixed  in  U,  this  stress  field  on  Y  shows  how  the  macroscopic  stress 
E(x)  =<  <r^(x,y)  >  is  localized  in  an  eY-period  (cell)  at  the  point  x  6  fl. 

This  procedxire  can  be  rigorously  justified  in  the  sense  that  it  can  be  shown  that 
when  €  tends  to  zero,  the  stress  field  a*{x)  tends  to  E(x).  Nevertheless,  a^{x,xf€) 
is  a  better  approximation  to  o*(x)  than  E(x).*  The  macroscopic  stress  field  E(x)  is 
just  a  mean  value  while  <7^(x,x/£)  takes  into  account  the  fine  periodic  structure  of 
the  composite  material. 


2.3  Comments 

These  results  demonstrate  the  effectiveness  of  homogenization  theory  in  comput¬ 
ing  the  mechanical  characteristics  of  composite  materials.  Validity  of  the  results 
depends  on  the  assumptions  made  on  the  shapes  and  layout  of  fibers  (especially  pe¬ 
riodicity).  Of  course,  a  method  which  permits  treatment  of  random  fiber  directions 
would  be  able  to  capture  the  actual  physical  situation  more  precisely.  However,  the 
undeniable  advantage  of  this  method  derives  &om  its  capability  to  supply  complete 
sets  of  parameter  values,  which  are  mutually  coherent.  Comparisons  with  exper¬ 
imental  evidence  and  other  techniques  such  as  the  Halpin-Tsai  formulas  indicate 
that  the  estimates  based  on  the  homogenization  theory  are  more  accurate  than 
those  based  on  other  methods  [7,15,16]. 

An  important  challenge  is  to  move  from  the  ability  to  characterize  the  material 
properties  of  composites  to  characterizing  the  macroscopic  behavior  of  structures 
fabricated  from  these  advanced  materials.  There  has  been  very  little  work  on  this 
subject,  and  it  is  important  that  the  software  modeb  provide  a  basb  for  the  devel¬ 
opment  of  CAD  toob  for  treatment  of  these  bsues.* 

'Precisely,  the  norm  convergence  implies  that  cr*(z)  —  <r^{x,  z/e)  tends  to  sero  for  almost 

every  point  in  0,  while  the  (weak)  £'(0)  convergence  governing  the  approach  of  o*  to  S  does  not. 

'See  [83]  for  a  general  discussion  of  the  impact  of  composites  on  design  of  high  performance 
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The  basic  analytical  method  can  be  extended  to  treat  elementary  structures  fab¬ 
ricated  from  composites.  In  the  next  section  we  discuss  the  treatment  of  laminates 
formed  from  fiber-reinforced  composites.  The  analysis  of  such  systems  has  been  a 
problem  in  the  past,  see  the  discussion  in  [71].  The  homogenization  method  can 
readily  treat  laminates  with  a  large  number  of  layers;  a  problem  that  has  presented 
difficulties  for  other  methods.  More  complex  structures  can  also  be  treated  by  it¬ 
eratively  applying  the  method.  For  example,  in  treating  tubular  members  made 
from  composites  one  can  compute  the  materiad  properties  of  the  composite  using 
the  homogenization  procedure,  evaluate  the  stiffness  and  bending  properties  of  the 
tube  based  on  this,  and  then  refine  the  model  by  examining  the  distribution  of 
microstresses  in  the  wall  of  the  tube  as  a  function  of  fiber  alignment  and  geometry. 
The  deformations  of  a  sleeve  coupling  element  may  be  similarly  treated. 

More  complex  molded  structures  are  a  different  matter,  since  the  macroscopic 
geometry  of  the  element  will  have  to  be  considered  even  in  the  local  analysis.  That 
is  the  “shape”  of  the  region  fl  will  be  very  complex  for  certain  molded  structures; 
and  it  will  have  to  be  treated  e3q)licitly  in  the  homogenization  analysb.  One  possible 
approach  is  to  use  an  iterated  finite  element  procedure  in  which  finite  elements  in 
the  large  scale  representation  of  the  complex  region  fl  are  themselves  represented  as 
“composite  elements”  which  are  are  in  ttim  treated  as  a  collection  of  “cells”  whose 
compliance  properties  are  evaluated  using  the  homogenization  procedure.  The  be¬ 
havior  of  the  overall  molded  structure  would  then  be  derived  by  a  “superposition” 
of  the  individual  behaviors  of  each  of  its  “composite  elements.”  Implementation  of 
this  procedme  would  be  a  very  challenging  problem  in  numerical  (finite  element) 
analysis.  The  simplicity  and  similarity  of  the  “composite  elements”  would  have  to 
be  exploited  in  setting  up  the  analysis.  If  properly  formulated,  the  analysis  of  the 
“composite  elements”  could  be  done  in  parallel. 

systems.  See  the  papers  [45,51,72]  for  discassions  of  design  objectives  for  laminated  fibrons  composite 
plates  and  related  structures  .  The  homogenisation  method  u  applied  to  design  in  [48]. 
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1  :  lut  period 
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by  0  double  fiber  rowt  forced  by  t  auUlple 
fiber  rowt. 


Figure  4:  Multiple  layer  reinforced  composites. 


3  Numerical  Evaluation  of  Effective  Moduli  and 
Microscopic  Stress  Fields 


3.1  Application  of  the  Method  to  a  Cross  Ply  Laminate 


To  illustrate  use  of  the  homogenization  method  in  a  concrete  case,  consider  a  cross 
ply  lamination  consisting  of  orthotropic  layers  with  an  alternating  orientation  of 
90"  between  the  layers.  We  assume  that  the  z  and  the  y  axes  are  parallel  to  the 
fiber  reinforcement  in  the  odd  and  the  even  numbered  layers  respectively.  All  the 
odd  layers  are  identical.  As  are  the  even  layers.  To  compare  the  results,  we  shall 
consider  three  types  of  layers:  (See  Figure  4.) 


Example  1: 
Example  2: 
Example  3: 


Each  layer  is  reinforced  by  a  single  fiber  row. 
Each  layer  is  reinforced  by  a  double  fiber  row. 
Each  layer  is  reinforced  by  multiple  fiber  rows. 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


SEI-TR-87-05 


Systems  Engineering,  Inc. 


17 


I 


Period  Y  (»1i>9''* 


Period  Y  (double  fiber  rew*) 


Figiire  5:  Triangular  finite  element  mesH  of  the  cross  section  of  a  period  Y. 


The  first  two  examples  are  solved  by  using  the  homogenization  procedure  in 
three  dimensions.  The  last  example  is  solved  in  two  steps  by  using  a  two-dimensional 
analysis: 


(i)  Homogenization  of  each  layer  reinforced  by  fibers  running  in  the  same  direction. 

The  problem  defined  by  (22)  is  bidimensional. 

(ii)  Homogenization  of  the  cross  ply  lamination  using  the  computed  results  for  each 

layer  given  by  the  first  step. 

The  munerical  results  are  obtained  by  using  a  finite  element  code.  See  Figure  5. 


3.2  Three  Dimensional  Homogenization  of  the  Cross  Ply 
Laminate 

To  obtain  the  homogenized  moduli  qff  (25),  it  is  necessary,  to  first  compute  the 
functions  which  are  solutions  of  elliptic  boimdary  value  problems  on  the  basic 
period  (22).  To  compute  the  functions  by  a  finite  element  code,  the  system  (22) 
is  reformulated  in  variational  form 

=  {t;  =  (t;i,V2,  Vs),®  e  {H*(y)}’,Si  is  y -periodic 
\  av(x*'ill)  =  fy  aijkt(y)ei,  (v)dj/,  Vv  e  V 


(30) 
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where 

The  last  integral  in  (30)  can  be  written  in  terms  of  the  surface  load: 

with 

•P? ^  =  [o-xie^nt  =  [onejn^ 

FP  =  [o3ja]n£  Fj”  =  [oajejni  F,**  =  [a22o]n£ 

Fi“  =  (assnjnr  F,”  =  (aasajnr  F^  =  [aasojnt 

Fi”  =  (ajanjne  F^  =  [a23£2]n£  Fj”  =  [onojn£ 

F/’  =  [aianjn/  F^^  =  (aiaajn/  F^^  =  [anojnf 

Fl^  =  [aiia\ni  F^  =  [oucjnr  F,”  =  [aijejn^ 

where  i%=  (ni,  nj,  ns)  is  the  outward  unit  normal  to  the  interface  of  the  components 
(fiber-resin)  and  where  the  bracket  {«]  denotes  the  jump  of  a  function  across  the 
interface. 

Given  numerical  representations  of  the  functions  ^he  homogenized  coef¬ 
ficients  qij-u  of  the  composite  system  are  the  mean  value,  over  a  period  Y  of  the 
corresponding  Oijkt  altered  by  a  corrector  term  (25)  depending  on  the  explicitly: 

f  qaiu  =<  a<,M(y)  >  -qf^u 
\  9i}U  = 

3.3  Two-Dimensional  Homogenization  of  the  Cross  Ply  Lam 
inate 

As  stated  previously,  the  Example  3  is  solved  in  two  steps: 

Step  1:  Homogenization  of  a  composite  reinforced  by  fibers  running  in  the  same 
direction.  See  Figiure  6. 

Calculation  of  the  homogenized  coefiScients  9,yju  requires  resolution  of  (22).  In 
the  present  case  the  coefficients  a,>M(y)  are  independent  of  y^.  Consequently,  the 
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Figure  6:  Composite  reinforced  by  fibers  running  in  the  same  direction 

fields  x^(y)  independent  of  ya.  In  (22)  (30)  those  indices  referring  to  dfdy^ 

give  zero  contribution.  This  makes  the  computation  of  x*'  (y)  ^  dimensional 
problem.  In  addition,  we  have: 

x“(y)  =  x“(yi.y3)  =  [xi*(yi,y2),x"(yi.y2),o] 

These  functions  are  therefore  solutions  of  a  plane  strain  elasticity  problem.  And, 

x"(y)  =  x"(yi,y2)  =  lo,o,x"(yi,y2)]  for  (fc,f)  =  [(i,3),(2,3)] 

These  two  functions  for  solutions  of  a  scalar  problem  in  (Details  of  this  reduction 
are  given  in  [7].) 

Step  £:  Homogenization  of  the  cross  ply  lamination.  See  Figure  7. 

Using  the  homogenized  moduli  of  each  layer  (from  step  one),  the  homogenization 
formulas  are  considerably  simplified.  The  problem  (22)  is  reduced  to  a  system  of 
differential  equations  which  may  be  solved  explicitly. 
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HiTtNPACE 


Figure  7:  Multiple  layers:  Each  layer  possesses  a  plane  of  elastic  symmetry  normal 
to  the  X3  axis. 

3.4  Sample  numerical  results 

3.4.1  Effective  Parameter  Models 

The  results  obtained  (from  a  finite  element  code)  for  the  two  first  examples  are 
compared  with  those  obtained  for  the  last  example  3  in  the  tables  which  follow: 
Each  component  is  assumed  elastic,  homogeneous,  and  isotropic  with: 

E/  =  84,000  MPa  Uj  =  ,22  for  the  fiber 
E,=  4,000  MPa  Ur  =  .ZA  for  the  resin 

Results  for  3  impregnations  of  resin  (36%,  50%,  65%)  are  shown  in  Tables  1,2,3. 


3.4.2  Numerical  Results  for  the  Microscopic  Stress-Field 

In  the  previous  paragraph,  we  have  shown  that  we  get  the  homogenized  moduli 
<liju  from  the  six  vector  fields  x*^(y)*  We  can  then  solve  the  homogenized  elastic 
problem  (27) (28),  on  0  which  gives  the  macroscopic  stress  field  £(x),  for  x  €  H. 
The  computations  of  the  microscopic  stress-field  and  stress  forces  at  the  interface 
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Example 

El 

(MPa) 

E2 

(MPa) 

E3 

(MPa) 

G23 

(MPa) 

G13 

(MPa) 

G12 

(MPa) 

1^23 

*^13 

1/12 

37,700 

37,700 

18,900 

5,600 

5,600 

6,775 

.27 

.27 

.135 

37,780 

37,780 

19,850 

5,570 

5,570 

6,670 

.259 

.259 

.137 

37,750 

37,750 

20,614 

5,298 

5,298 

6,304 

.24 

.24 

.138 

Table  1:  Comparative  table  for  36%  resin  by  volume. 


Example 

El 

(MPa) 

E2 

(MPa) 

E3 

(MPa) 

G23 

(MPa) 

G13 

(MPa) 

G12 

(MPa) 

U23 

1/13 

1/12 

1 

28,340 

28,340 

12,996 

3,800 

3,800 

4,300 

.32 

.32 

.122 

2 

28,400 

28,400 

13,280 

3,780 

3,780 

4,276 

.308 

.308 

.124 

3 

28,801 

28,801 

13,563 

3,598 

3,598 

4,168 

.297 

.297 

.124 

Table  2:  Comparative  table  for  50%  resin  by  volume. 


Example 

El 

(MPa) 

E2 

(MPa) 

E3 

(MPa) 

G23 

(MPa) 

G13 

(MPa) 

G12 

(MPa) 

1/23 

1/13 

1/12 

1 

20,500 

20,500 

9,471 

2,880 

2,880 

3,076 

.33 

.35 

.126 

2 

20,560 

20,560 

9,565 

2,850 

2,850 

3,057 

.35 

.35 

.128 

3 

20,237 

20,237 

9,243 

2,657 

2,657 

2,910 

.35 

.35 

.126 

Table  3:  Comparative  table  for  65%  resin  by  volume. 
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between  fiber  and  resin,  are  particularly  important  because  they  can  initiate  cracks 
and  delaminations.  The  localization  procedure  allows  an  easy  computation  of  these 
microscopic  stress-field  and  stress  forces. 

Numerical  results  for  an  unidirectionally  fiber-reinforced  composite  subjected  to 
a  shearing  stress-field  normal  to  the  direction  of  the  fibers  are  shown  in  Figure  8. 
Other  cases  can  be  treated  in  a  similar  manner.  The  stress  forces  a  *  n  at  the 
interface  and  the  components  of  microscopic  stress  field  are  plotted.  The  density  of 
lines  plotted  indicate  the  level  of  microscopic  stress  concentration. 

The  elastic  moduli  of  the  components  are: 

Fiber:  E  =  84000  Mpa  v  =  0.22 

Resin:  E  =  4000  Mpa  v  =  0.34 


The  corresponding  homogenized  moduli  for  50%  resin  impregnations  are 


El  =  1014lMpa 

t/Z2 

=  0.287 

G32  —  3106 

E2  =  96S5Mpa 

Uzi 

=  0.281 

G31  =  3386 

EZ  =  35655Afpa 

Ui2 

=  0.353 

Gi2  =  2606 

4  Numerical  Results  for  Eflfective  Moduli  and  Mi¬ 
croscopic  Stress  Fields 


In  this  section  we  present  transcripts  of  interactive  sessions  with  the  NeMCAP  software 
system.  Computation  of  the  effective  moduli  is  treated  in  the  next  subsection. 
Evaluation  of  the  microscopic  distribution  of  stresses  at  the  fiber-matrix  interface 
is  treated  in  the  last  subsection.  In  Listings  1  and  2  at  the  end  of  the  report  we 
provide  more  extensive  examples  of  the  performance  of  the  program  modules. 


4.1  EjSective  Moduli  for  Various  Fiber-Matrix  Geometries 

In  the  session  which  follows  the  homogenized  compliance  tensor  defining  the  re¬ 
lationship  between  the  stress  and  strain  tensors,  {e}  =  [E]  {o},  is  computed  for 
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sample  fiber  types  and  orientations.  Specifically,  we  computed  the  system 


Cll 

r  1 

El 

El 

0 

0 

0  ■ 

On 

^22 

—  Hix 

El 

1 

0 

0 

0 

s 

022 

*33 

El 

_iauL 

i  0 

0 

0 

O33 

«23 

0 

0 

0  2*1 

0 

0 

023 

ei3 

0 

0 

0 

0 

1 

2Git 

0 

Oi3 

*12 

0 

0 

0 

0 

0 

1 

iGit  j 

<Ti2 

Here  =  1,2,3  are  Young’s  moduli  for  the  fiber  (composite)  directions  x,-,i  = 
1,2,3;  G^,-,*,y  =  1,2,3  are  the  fiber  (composite)  shear  moduli,  and  =  1,2,3 

are  the  fiber  (composite)  Poisson  coefficients. 

The  homogenization  methodology  was  used  to  compute  the  effective  parame¬ 
ters.  As  shown  in  the  previous  section,  this  requires  solution  of  an  elliptic  partial 
differential  equation  in  a  typical  "cell”  of  the  (periodic)  structure.  This  solution 
was  found  by  a  finite  element  procedure.  The  solution  is  then  used  in  an  averaging 
procedure  which  produces  the  effective  parameters. 


4.1.1  Design  of  the  MeMCAP  Supervisor 

We  have  developed  a  menu-driven  supervisor  program  called  the  "Metal-Matrix 
Composites  Analysis  Program  (NeNCAP)”  to  facillitate  interaction  with  the  nximer- 
ical  programs  which  compute  the  effective  parameters  of  the  "homogenized”  repre¬ 
sentation  of  fiber  reinforced  materials.*^  An  overview  of  this  software  is  shown  in 
Figure  9. 

The  fiow  diagram  of  the  menu-oriented  supervisor  module  is  shown  in  Figure  10. 

Currently,  the  software  allows  the  user  to  (i)  get  help,  (ii)  set  up  a  problem  via 
a  simple  screen-oriented  editor  and  (iii)  submit  a  batch  job  to  compute  the  effective 
parameters  of  the  current  problem.  The  numerical  results  are  saved  in  a  .LOG  file. 

fact,  most  of  the  data  used  in  the  computations  here  is  for  resin  based  composites.  We  have 
not  used  metal  matrix  data  (in  this  unrestricted  document)  in  part  to  avoid  issues  related  to  the 
ITAR. 
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ran 


Figure  0:  Overview  of  MeMCAP  software  system. 
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Figure  10:  Flow  diagram  of  MeMCAP  Supervisor 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


SEI-TR-87-05  Systems  Engineering.  Inc.  27 

The  system  is  written  in  Fortran.  It  assumes  the  user  has  a  VTlOO  terminal  or 
compatible.  The  system  is  just  a  prototype;  we  plan  significant  enhancements  in 
further  work. 


4.1.2  Sample  Session  with  the  MeMCAP  Interface 


The  output  horn  a  sample  run  is  shown  below.  Each  page  corresponds  to  a  display 
screen  on  the  user’s  terzmnal. 


$  run  mameap 


WELCOME 

TO 

*tm****************************** 

*  MaNCAP  * 

*  Natal  Matrix  Conpoaitas  * 

*  Analysis  Program  * 

a******************************** 


PROTOTYPE  PROGRAM  FOR  COMPUTING  THE  EFFECTIVE  PARAMETERS 
OF  METAL  MATRIX  COMPOSITES 

DEVELOPED  BY 

SYSTEMS  ENGINEERING  INCORPORATED 
7833  WALKER  DRIVE.  SUITE  308 
GREENBELT.  HD  20770 


PRESS  THE  "RETURN"  KEY  TO  CONTINUE  ♦** 
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MAIN  MENU 


ABBREVIATIONS 


1  HELP  (BRIEF  DESCRIPTION  OF  SOME  ITEMS)  -  (HELP) 

2  PREPARE  DATA  FILE  FOR  COMPUTATIONS  —  (PREP) 

3  COMPUTE  EFFECTIVE  PARAMETERS  -*  (COMP) 

4  . 

6  . 

6  . 

0  QUIT 

PLEASE  ENTER  TOUR  SELECTION  NUMBER: 
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HELP  -  Help  available  on  the  following  topics: 

HELP  -  displays  this  list.  HISTORY  -  history  of  this  project. 

NEWS  -  info,  on  revisions.  COMP  -  confute  elf.  paraneters 

PREP  -  preparation  of  data  file. 

SUMMARY  '  general  description. 


1.  TO  SEE  A  DIRECTORY  OF  HELP  ITEMS.  TYPE  "HELP". 

2.  OR  ENTER  NAME  OF  ITEM  FOR  WHICH  YOU  WANT  HELP. 

3.  WHEN  FINISHED.  TYPE  "EXIT"  OR  "QUIT". 


EXIT 
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MAIN  MENU 


ABBREVIATIONS 


1  HELP  (BRIEF  DESCRIPTION  OF  SOME  ITEMS)  -  (HELP) 

2  PREPARE  DATA  FILE  FOR  COMPUTATIONS  —  (PREP) 

3  COMPUTE  EFFECTIVE  PARAMETERS  —  (COMP) 

4  . 

5  . 

6  . 

0  QUIT 

PLEASE  ENTER  TOUR  SELECTION  NUMBER: 
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ITEM  DESCRIPTION 

1.  Project  title 

2.  Project/Contract  Number 

3.  Generic  Name 

4.  IFIB 


TYP  CURRENT  VALUE 

C  Metal  Matrix  Compoaites  Analysis 
C  test 
C  TEST 
1  1 


Fiber  type: 

1  »  Isotropic  Circular  Fiber 

2  «  Orthotropic.  Circular  Fiber 

3  s  Isotropic  Kidney  Fiber 

4  -  Orthotropic  Kidney  Fiber 

5  s  Isotropic  Staggered  Fiber 


TTP:  C^CHARACTER,  I*INTEGER.  L«L0GICAL.  R»REAL.  XeCOMPLEX 


1.  TO  CHANGE  CURRENT  ITEM.  ENTER  NEW  VALUE.  THEN  BIT  RETURN. 

2.  OR  ENTER  *«"  FOLLOWED  BY  ITEM  NUMBER  YOU  WANT  TO  CHANGE. 

3.  TO  WRITE  BACK  FILE  AND  EXIT.  TYPE  "EXIT".  ELSE  "QUIT". 


31 


QUIT 
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MAIN  MENU 


ABBREVIATIONS 


1  HELP  (BRIEF  DESCRIPTION  OF  SOME  ITEMS)  -  (HELP) 

2  PREPARE  DATA  FILE  FOR  COMPUTATIONS  —  (PREP) 

3  COMPUTE  EFFECTIVE  PARAMETERS  —  (COMP) 

4  . 

6  . . 

6  . 

0  QUIT 

PLEASE  ENTER  YOUR  SELECTION  NUMBER: 

MAIN>3 

SUBMITTING  BATCH  JOB.  LOGFILE  IS  ••  COMACH.LOG" 

Job  COMACH  (queue  SYS$BATCH.  entry  1098)  started  on  SYS$BATCH 


The  results  for  this  case  are  shown  below  (i.e.,  the  following  is  a  listing  of  the 
.LOG  file  created  by  the  module  COMACH.EXE): 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


SEI-TR-87-05 


Systems  Engineering.  Inc. 


33 


$RUN  SEI$USR: [CRANE] CONACH 
GENERIC  NAME  : 

IFIB  ?  (  «  1  ISOTROPIC  CIRCULAR  FIBER) 

(  *  2  ORTHOTROPIC  CIRCULAR  FIBER) 

(  «  3  ISOTROPIC  KIDNEY  FIBER) 

(  «  4  ORTHOTROPIC  KIDNEY  FIBER) 

(  »  5  ISOTROPIC  STAGGERED  FIBER) 

(  «  6  ORTHOTROPIC  STAGGERED  FIBER) 

(  *  7  HEXAGONAL  CELL  ISOTROPIC  KIDNEY  FIBER) 

(  s  8  HEXAGONAL  CELL  ORTHOTROPIC  KIDNEY  FIBER) 

1 

lOPT  *  ?  IF  »  1  PROVIDE  THE  SIDE  RATIO  (  SIDE_//_Y  /  SIDE_//_X  ) 

AND  THE  RESIN  RATIO 

IF  s  2  PROVIDE  THE  SIDE  RATIO  .  THE  RESIN  RATIO, 

HDIST  (  MESH  DENSITY  PARAMETER  ) 

1 

SIDE  RATIO  ?  (REAL) 

1.0000000 

RESIN  RATIO  ?  (  >  0.2348894  ) 

O.BOOOOOO 

E  .  NU  «  ?  (YOUNG  MODULUS  POISSON  COEFF.  ISOTROPIC  FIBER) 
84000.0000000 
0.2200000 

E  .  NU  «  ?  (YOUNG  MODULUS  POISSON  COEFF.  ISOTROPIC  RESIN) 
4000.0000000 
0.3400000 
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THE  FIBERS  ARE  PARALLEL  AT  XI 

**********4i4>*4t*******4i**4>4<*4<4i*** 


NOTATIONS 


E1.E2.E3:  YOUNG  MODULI  FIBER  (COMPOSITE)  DIRECTIONS  X1.X2,X3 

GIJ  :  FIBER  (  COMPOSITE  )  SHEAR  MODULI 

NU(IJ)  :  FIBER  (  COMPOSITE  )  POISSON  COEFFICIENT 

ER  :  RESIN  YOUNG  MODULUS 

NUR  :  RESIN  POISSON  COEFFICIENT 

NEL  :  NUMBER  OF  ELEMENTS 

NOE  :  NUMBER  OF  NODES 

TXR  :  RESIN  IMPREGNATED  RATIO  IN  VOLUME 


****************^t******^^**********^l^**%^^:lt:*^L1^*^t****************************l^ 

*  ISOTROPIC  CIRCULAR  FIBER  ♦ISOTROPIC  RESIN^  MESH  ♦ 


•  El  =  84000.  .  612  «  34426.  .  NU12  =  0.220  ♦  ER  *  4000.  *  NEL  =  176  ♦ 

•  E2  =  84000.  .  G13  =  34426.  .  NU13  =  0.220  ♦  ♦  NOE  =  106  ♦ 

•  E3  =  84000.  .  G23  =  34426.  .  NU23  «  0.220  ♦  NUR  »  0.340  ♦  TXR  =  .5082  ♦ 


*******m*******************0***0*****nmt********************************4L}tm** 
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*il,tHHr*1H,m***********m******************************************************* 

*  HOMOGENIZED  ELASTIC  TENSOR  * 

*  . • 


♦ 

1  1 

2  2 

3  3 

2  3 

1  3 

1  2 

* 

0.462BE-P0B 

0.B296E+04 

0.B297E404 

-0.7273E-13 

O.OOOOE+00 

O.OOOOE+00 

« 

0.B296E-t-04 

0.1468E-f0B 

0.4762E+04 

-0.3214E-12 

O.OOOOE+00 

O.OOOOE+00 

« 

0,B297E+04 

0.4762E+04 

0.1468E+0B 

-0.8976E-13 

O.OOOOE+00 

O.OOOOE+00 

-0.7273E-13 

-0.3214E-12 

-0.8976E-13 

0.3113E-f04 

O.OOOOE+00 

O.OOOOE+00 

* 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE-pOO 

0.4077E+04 

0.9916E-04 

* 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.9916E-04 

0.4077E+04 

**************i^***^*,tt^t***m**************************m********************** 


***4>***4i*****4<**4<4<4t****4i4>4<<*4i4<4>****4>**4>4>*****************4>********4i*«>*4t4i4t*4>4i* 


* 

*  . 
* 
* 
* 
* 
* 
* 
* 


HOMOGENIZED  COMPLIANCE  TENSOR 


* 

* 

* 

* 

* 

* 

* 

* 


11  22  33  23 

0.2306E-04  -0.6283E-05  -0.6283E-05  -0.2010E-21 
-0.6283E-0B  0.7786E-04  -0.2299E-04  0.7228E-20 
-0.6283E-05  -0.2299E-04  0.778BE-04  -0.27B4E-21 
-0.2910E-21  0.7228E-20  -0.27B4E-21  0.3212E-03 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 


1  3 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 


1  2 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 


0.24B3E-03  -0.B967E-11 
•O.B967E-11  0.24B3E-03 


El 

0.4336B9D-P0B 

E2 

0.128441D-P0B 

E3  > 

0.12844BD-i-0B 

NU23  « 

0.29B2920-P00 

NU12  « 

0.272460D+00 

NU13  « 

0.2724B8D+00 

G23  « 

0.311300D+04 

G12  > 

0.4076B7D+04 

G13  = 

0.407662D-E04 

*************************’****L**m*i¥*******0*************iHlf******************* 
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ROTATION  45  DEGREES  AROUND  II 

*4i******4i***«***«*4i*«4t4>*«4i****4i 


****4>*4i4c*«**4i**4>****«******«******«******4i********4i*«4i4t****4i*****«**4i**4>**4i* 


*  HOMOGENIZED  COMPLIANCE  TENSOR  (BISECTING  DIRECTIONS)  * 

•  . . . * 

*11  22  33  23  13  12  * 


*  0.2306E-04  -0.6283E-0B  -0.6283E-05  0.6860E-10  O.OOOOE+00  O.OOOOE+OO  * 

*  -0.6283E-05  0.1077E-03  -0.5288E-04  -0.1208E-08  O.OOOOE+OO  O.OOOOE+OO  ♦ 

*  -0.6283E-05  -0.B288E-04  0.1077E-03  -0.1208E-08  O.OOOOE+OO  O.OOOOE+OO  ♦ 

*  0.6650E-10  -0.1208E-08  -0.1208E-08  0.2017E-03  O.OOOOE+OO  O.OOOOE+OO  • 

*  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.2453E-03  -0.1B74E-08  ♦ 

*  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  -0.1574E-08  0.2453E-03  ♦ 

********««***4<**«4i4>*****«*******4i*«****«***«***«*4t***4t4>»4i****4i************** 


* 

*  El  «  0.4336B90+05 

* 

*  NU23  *  0.490768D+00 

* 

*  G23  -  0.49B8050+04 

* 


E2  «  0.928151D+04 
NU12  s  0.2724B9D+00 
G12  «  0.4076B9D+04 


* 

E3  e  0.9281B1D+04  * 

* 

NU13  =  0.272459D+00  * 

* 

G13  «  0.4076B9D+04  * 

* 
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FORTRAN  STOP 

CRANE  Job  terminated  at  ll-MAY-1087  14:34:58.43 


Accounting  Inlormatlon: 


Buffered  1/0  count: 

08 

Peak  working  set 

size: 

1354 

Direct  1/0  count: 

164 

Peak  page  file  size: 

3151 

Page  faults: 

1363 

Mounted  volumes: 

0 

Charged  CPU  time: 

0  00:01:07.40 

Elapsed  time: 

0  00 

.‘01:20.34 

The  program  was  run  on  SEI’s  VAX  11/750,  which  has  a  floating  point  accel¬ 
erator  and  6  Megbytes  of  core  memory.  Listing  1  at  the  end  of  the  report  contains 
analyses  of  eight  different  fiber-matrix  geometries  to  illustrate  the  range  of  options 
available  in  this  module  of  the  system. 
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4.2  Evaluation  of  the  Microscopic  Stress  Fields 

The  MeMCAP  software  includes  a  module  (C0MPQ.EXE)  for  the  evaluation  of  micro¬ 
scopic  stress  distributions  in  a  typical  cell  of  the  composite  material.  The  system 
is  setup  to  handle  cases  in  which: 

1.  The  fibers  are  run  in  the  same  direction  and  are  periodically  distributed  with 
a  (i)  staggered  or  (ii)  aligned  rectangular  base  period. 

2.  The  fiber  material  is  (i)  isotropic  or  (ii)  orthrotropic. 

3.  The  matrix  material  is  (i)  isotropic,  (ii)  orthrotropic,  or  (iii)  incompressible. 

4.  The  fibers  are  (i)  circular  or  (ii)  elliptical.'^ 

5.  The  microscopic  stress  field  is  computed  for  six  main  (stress)  loadings:  (i) 
Sn  =  hSii  =  0,t  5^  i,  (ii)  $22  =  1,5,-,.  =  0,»  ^  i,  (iii)  533  =  hSii  =  0,*V  y, 
(iv)  523  =  1,  Si,-  =  0, 1 y,  (v)  5i3  =  1,  Sij  =  0,  *  ,#  y,  (vi)  5i2  =  1, 5y  =  0,  t  # 

y* 

The  finite  element  mesh  is  automatically  generated  given  the  description  of 
the  cell  geometry  and  the  fiber-matrix  impregnation  ratio.'*  A  ratio  is 

provided  for  each  fiber  shape  and  geometry. 

The  software  requires  the  following  input  data: 

1.  Output  device:  terminal  or  filename. 

2.  Fiber  characteristics:  aligned  or  staggered;  isotropic  or  orthotropic;  spacing 
of  fibers  in  x  and  in  y;  and  shape  -  circular  or  ellipse  (a,&,n}. 

3.  Fiber-matrix  impregnation  ratio. 

4.  Elasticity  characteristics:  Matrix  -  isotropic  (E,i/);  orthotropic  (Ei,Gi,-,i/,-,); 
or  incompressible  {E,v);  and  Fiber  -  isotropic  (E,  u)  or  orthotropic  (E,- ,  G,-,- ,  i/,y) . 

Actually,  a  generaliaed  elliptical  shape  is  permitted  (x/o)"  -h  (y/h)"  =  1. 

optional  parameter  (HOIST)  can  be  used  to  refine  the  mesh  near  the  boundary. 
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The  software  produces  the  following  output  information  either  on  the  user’s 
terminal  or  into  a  file,  as  directed  by  the  user: 

1.  A  summary  of  the  input  data. 

2.  The  homogenized  elasticity  tensor. 

3.  The  homogenized  compliance  tensor. 

4.  Equivalent  elastic  moduli. 

5.  A  stimmary  of  the  finite  element  mesh  including  number  of  elements  and 
number  of  nodes.  For  each  element  the  program  provides:  number,  material 
number  (1  for  fiber,  2  for  matrix),  node  number,  and  node  coordinates.  This 
informationn  serves  to  localize  the  microstresses  in  the  fiber-matrix  system. 

6.  For  each  stress  loading  the  microscopic  stress  field  by  element  ntimber  (in  the 
mesh). 

7.  The  stress  force  at  the  interface  between  znatrix  and  fiber  in  terms  of  element 
number,  node  coordinates,  and  stress  force  components. 

The  results  of  a  typical  interactive  session  with  the  software  are  given  in  Listing 
2.^^  The  listing  of  microscopic  stress  distributions  by  mesh  element  is  tiseful  in 
anticipating  design  objectives  and  limitations.  It  is  also  possible  to  have  the  results 
displayed  in  a  graphical  format.  Figxires  11,12,13  show  the  microstress  distribution 
in  a  typical  cell  in  response  to  various  loadings.  The  colors  correspond  to  various 
stress  levels. 

^*To  save  space,  we  omit  details  of  the  interaction  with  the  system  through  the  supervisor  interface. 
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Figure  11:  Microstress  distribution  in  a  typical  cell. 
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Figure  13:  Microstress  distribution  in  a  typical  cell  (alternate  geometry). 
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5  Conclusions  and  Further  Research  and  Devel¬ 
opment 

5.1  Accomplishments  of  the  Project 

The  objective  of  oiur  Phase  I  research  program  was  to  demonstrate  the  feasibility 
of  the  homogenization  method  to  compute  the  effective  parameter  representations 
of  fiber  reinforced  composites.  We  have  produced  analytical  results  and  a  software 
system  which  meet  this  objective.  It  is  important  to  stress  that  the  analytical 
methods  are  systematic  and  consistent. 


5.2  Further  Phase  II  Research  and  Development 

Several  enhancements  to  the  software  are  possible  which  would  male  it  a  more  com¬ 
plete  CAD  system  for  investigation  and  design  of  composite  materials  and  structures 
constructed  from  such  materials. 


1.  As  we  have  already  indicated,  it  would  be  relatively  straight-forward  to  add 
enhancements  to  treat  the  viscoelastic  and  thermal  properties  of  (periodic) 
composites  (in  the  linear  regime). 

2.  Similarly,  it  is  possible  to  adapt  the  methods  to  treat  plasticity  in  composites 
-  a  particularly  important  consideration  in  the  analysis  and  design  of  metal 
matrix  composites. 

3.  It  is  possible  to  adapt  the  analytical  methods  to  treat  both  short-fiber  and 
particulate  reinforced  composites;  however,  in  these  cases  it  is  important  to 
be  able  to  treat  systems  with  a  random  infrastructure. 

4.  Systems  with  a  random  infrastructure  can  be  treated  using  homogenization. 
The  general  theory  has  been  worked  out  [59].  Extensive  applications  have 
been  given  both  for  nonlinear  problems  [10]  and  for  systems  with  a  discrete 
infrastructure  (heat  propagation  on  a  random  lattice)  [50].  In  related  research 
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SEI  is  developing  numerical  methods  for  the  treatment  of  homogenization 
problems  for  random  structures.^^ 
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Listing  1:  Computation  of  Effective  Moduli 


|lUN  COHACH  > 

l^£NERZC  NAME  : 

?FIB  ?  <  =  1  ISOTROPIC  CIRCULAR  FIBER) 

i  =  2  ORTHOTROPIC  CIRCULAR  FIBER) 

<  =  3  ISOTROPIC  KIDNEY  FIBER) 

C  =  A  ORTHOTROPIC  KIDNEY  FIBER) 

C  =  5  ISOTROPIC  STAGGERED  FIBER) 

(  =  6  ORTHOTROPIC  STAGGERED  FIBER) 

<  =  7  HEXAGONAL  CELL  ISOTROPIC  KIDNEY  FIBER) 

(  s  S  HEXAGONAL  CELL  ORTHOTROPIC  KIDNEY  FIBER) 

loPT  s  ?  IF  =  1  PROVIDE  THE  SIDE  RATIO  (  SIDE_//_Y  /  SIDE_//_X  ) 

AND  THE  RESIN  RATIO 

IF  =  2  PROVIDE  THE  SIDE  RATIO  #  THE  RESIN  RATIO^ 

HDIST  (  P.ESH  DENSITY  PARAMETER  ) 

1 

SIDE  RATIO  ?  (REAL) 

1 

RESIN  RATIO  ?  (  >  0.2348394  ) 

-5 

E  NU  -  ?  (YOUNG  MODULUS  POISSON  COEFF.  ISOTROPIC  FIBER) 

84000.  .22 

£  #  NU  »  ?  (YOUNG  MODULUS  POISSON  COEFF.  ISOTROPIC  RESIN) 

4000.  .34 


THE  FIBERS  ARE  PARALLEL  AT  XI 
********  ********  ************** 


NOTATIONS 

********* 


E1*E2#E3:  YOUNG  MODULI  FIBER  (COMPOSITE)  DIRECTIONS  X1*X2*X 

GIJ  :  FIBER  (  COMPOSITE  )  SHEAR  MODULI 

NL(IJ)  ;  FIBER  (  COMPOSITE  )  POISSON  COEFFICIENT 

cR  :  RESIN  YOUNG  MODULUS 

NUR  :  RESIN  POISSON  COEFFICIENT 

NEL  :  NUMBER  OF  ELEMENTS 

NOE  :  NUMBER  OF  NODES 

TXR  :  RESIN  IMPREGNATED  RATIO  IN  VOLUME 


* 

* , 

ISOTROPIC 

CIRCULAR 

FIBER 

♦ISOTROPIC 

_ _ 

RESIN« 

MESH 

* 

El  = 

84000.  . 

G12  = 

34426.  . 

NU12  = 

0.220 

« 

£R  s 

4000.  * 

NEL  > 

176 

* 

£2  = 

84000.  . 

G13  = 

34426.  . 

NU13  = 

0.220 

* 

* 

NOE  = 

105 

* 

E3  = 

84000.  . 

623  = 

34426.  . 

NU23  = 

0.220 

* 

NUR  = 

0.340  * 

TXR  = 

.5082 

*************************************************************************** 

*************************************************************************** 

*  HOMOGENIZED  ELASTIC  TENSOR  ^  j 

‘ . . 2“2 . 3"  3 . 2“  3 . rs . 1*2 . 


.  REPRODUCED  AT  GOVERNMENT  EXPENSE  ««««««««« 


0-4625E+05  0-52V6E+04  0.5297E+04  -0.7273E-13  O-COOOE^OO  O.OOOOE^OO 


Q.5296E+04  O-lAbSE-t-OS  0.4762E  +  Q4  -Q.3214E-12 
0.5297E+04  0-4762E+04  0.1468E+05  -0.8976E-13 


Q.Q000E400 

O.OOOOE'i-00 


Q.QQ00E400 

O.gOOOE^OO 


*  -0.7273E-13  -0-3214E-12  -0.8976E-13  0.3113E*04  O-OOOOE+00  O.OOOOE+OO 

*  O.OOOOE't-00  O.OOCOEi-OO  O.COOOE^OO  O.OOOOE4-00  0.4077E«^04  0.9916E-04 

*  O-OOOOE+00  O.OOOOE+00  O.COOOE+00  O.OOOOE+00  0.9916E-04  0-4077E+04 

*************************************************************************** 


*************************************************************************** 

*  HOMOGENIZED  COMPLIANCE  TENSOR 

*  . 2“! . 3*3 . 2*3 . 1*3 .  1  2 

*  0.2306E-04  -0-6283E-05  -0.6283E-05  -0.2910E-21  O.COOOE+OQ  O.OOOOE-fOO 

*  -0-6283E-05  0-7786E-04  -0.2299E-04  0.7228E-20  O.COOOE  +  00  O.OOOOE-EOO 

*  “0.6283E-05  “0-2299E-04  0-7785E-04  -0.2754E-21  O-OOOOE+00  O.OOOOE+OO 

*  -0-2910E-21  0-7228E-20  -0.2754E-21  0.3212E-03  O.OOOOE+00  O.OOOOE+00 

*  O.OOOOE+00  O.OOOaE+OO  O.COOOE+OO  O.OOOOE+00  0.2453E-C3  -0.5967E-11 

*  O-OOOOE+00  O-OOOOE+OO  O-OOOOE+00  O.OOOOE+00  -0.5967E-11  0.2453E-03 

*************************************************************************** 

**************************************************************************** 
«  * 

*  El  =  0.433659D+05  E2  =  0-12S4410+05  E3  =  C.  1 28445D't-05  * 

*  * 

*  NU23  =  0.295292D+C0  NL‘12  =  0.272460D+00  NUI 3  =  C-272458O+00  * 

*  * 

*  G23  =  0.31130QD+C4  612  =  0-4076570+04  613  =  0-407662D+04  * 

*  * 

************************  ************************  ************  It*  ********** 


ROTATION  45  DEGREES  AROUND  X1 

******************************* 


****************** 

*  HOKCGENl 


********************************************************** 
ZED  COMPLIANCE  TENSOR  (BISECTING  DIRECTIONS)  * 


*  11  2  2 

*  0-2306E-04  -0.6283E-05 

*  -0.6283E-05  0.1 077E-03 

*  -0.6283E-05  -0.52b8E-04 

*  0.6650E-10  ”0.1 203E-03 

*  O.OOOOE+00  O.OOOOE+00 

*  O-OOOOt+OO  O.OOOQ'+OO 
************************** 


3  3 

-0-6283E”05 


2  3 

0.6650E-1  0 


-0.3288E-04  -0-1208E-03 
C.1077E-03  -0.1208E-08 


-0.1203E-03 
O.QOOOE+00 
0. 00006+00 


0.2017E-03 

O.OOOOE+00 

O.COOOE+00 


1  3 

,00DOE  +  00  C. 
,C000E  +  00  0. 

,C000E+00  0. 

.OOOOE+00  0. 
,24536-03  ”0. 
,1574E-08  0. 


1  2 

OOOOE+00 

OOOOE+00 

OOOOE+00 

OOOOE+00 

1574E-08 

2453E-03 


************************** 

* 

*  El  =  0.4336590+05 

* 

*  NU23  =  C-490768D+C0 

* 

*  623  =  0.495605D+C4 

* 

************************** 
GENERIC  NAME  : 


************************************************** 

E2  =  0-928151D+J4  E3  =  C.928151D+04  * 

NU12  =  0-272459D+00  NU1 3  =  Q-272459D+00  * 

312  =  0-4076590+04  613  =  C.407659D+04  * 

* 

************************************************** 


IFIB  ?  (  -  1  ISOTROPIC  CIRCULAR  FIBER) 

(  s  2  ORTHOTROPIC  CIRCULAR  FIBER) 

(  s  3  ISOTROPIC  KIDNEY  FIBER) 

<  =  4  ORTHOTROPIC  KIDNEY  FIBER) 

(  =  5  ISOTROPIC  STAGGERED  FIBER) 

<  =  6  ORTHOTROPIC  STAGGERED  FIBER) 

(  =  7  HEXAGONAL  CELL  ISOTROPIC  KIDNEY  FIBER) 

C  =  8  HEXAGONAL  CELL  ORTHOTROPIC  KIDNEY  FIEER) 

lOPT  =  ?  IF  =  1  PROVIDE  THE  SIDE  RATIO  (  SIDE.//_Y  /  SIOE.//_X  ) 

AND  THE  RESIN  RATIO 

IF  =  2  PROVIDE  THE  SIDE  RATIO  ^  THE  RESIN  RATIO^ 

HDIST  (  P.ESH  DENSITY  PARAMETER  ) 

1 

SIDE  RATIO  ?  CREAL) 

1 

RESIN  RATIO  ?  <  >  0.2348894  ) 

-5 

El  ^  £2  ^  E3  =  ?  CORTHOTROPIC  FIBERS  //  A  XI) 

380000.  14300.  14300. 

612  ^  613  #  623  =  ?  CORTHOTROPIC  FIBERS  //  A  XI) 

38000.  38000.  20000. 

NU12  ^  NU13  /-  NU23  =  ?  CORTHOTROPIC  FIBERS  //  A  XI) 

.22  .22  .23 

E  #  NU  =  ?  CY0JN6  MODULUS  POISSON  COEFF.  ISOTROPIC  RESIN) 

3320.  .38 


THE  FIBERS  ARE  PARALLEL  AT  XI 
******************************** 


NOTATIONS  £1^E2/E3:  YCUNG  MODULI  FIBER  CCOMPOSITE)  DIRECTIONS  X1#X2# 

*********  oIJ  :  fiber  C  COMPOSITE  )  SHEAR  MODULI 

NUCIJ)  ;  FIBER  C  COMPOSITE  )  POISSON  COEFFICIENT 

ER  :  RESIN  YOUNG  MODULUS 

NUR  ;  KESIN  POISSON  COEFFICIENT 

NEL  :  number  OF  ELEMENTS 

NOE  :  NUMBER  OF  NODES 

TXR  :  RESIN  IMPREGNATED  RATIO  IN  VOLUME 


************  ******  ********************************************************* 

*  ORTHOTROPIC  CIRCULAR  FIBER  *ISOTROP1C  RESIN*  MESH 


*  El  =  38C000-  -  612  =  38000.  .  NU12  =  0.220  *  ER  *  3520-  *  NEL  =  176 

*  E2  =  14500.  .  613  =  33000.  -  NU13  =  0.220  *  *  NOE  =  105 

*  E3  =  14500-  -  623  =  20000-  -  NL'23  =  0.250  *  NUR  *  0.360  *  TXR  =  -5082 

*************************************************************************** 

********************************************  ****  **************  ************ 

*  HOMOGENIZED  ELASTIC  TENSOR 


REPRODUCED  AT  GOVERNMENT  EXPENSE  **♦♦*«*«« 


»  » 


* 

0-1912E+06 

0-4126E+04 

0.4126E+04 

-0-2165E-1 4 

O-OOOOE+OO 

O.OOOOE+00 

* 

0-4126E+04 

0-96b3E+04 

0.4199E+04 

0-3  317E-1  3 

O.OOOOE+00 

O.OOOOE+00 

* 

0-4126E+04 

0-4199E+04 

0-9684E+04 

-0-2809E-1  3 

O-COOOE+00 

O-OOOOE+OO 

* 

-0-2165E-14 

0-3317E-13 

~0-2809E“13 

0-2600E+04 

O.COOOE+00 

O-OOOOE+OO 

* 

O-OOOOE+00 

O-OOOOE+00 

O-OOOQE+OO 

O-OOOOE+OO 

0.3S70E+04 

0.1208E-02 

* 

O.OOGOE+00 

O-OOOOE+OO 

O-OOOOE+CO 

O-OOOCE+00 

0-120BE-02 

0-3570E+04 

*************************************************************************** 


************************* 

HOMOG 


************************************************* 

ENIZED  COMPLIANCE  lENSCR 


* 

1  1 

2  2 

3  3 

2  3 

1  3 

1  2 

* 

0-5299E-05 

-0-1  575E-05 

-0-1575E-05 

0.7<»95E-23 

O-OOOOE+OO 

O-OOOOE+OO 

* 

-0-1575E-05 

0-1 277E-03 

-0.5469E-04 

-0.2221E-20 

O-COOOE+00 

O-OOOOE+OO 

* 

-0-1 575E-05 

-G.5469E-04 

0.1277E-03 

0.2076E-20 

O-COOOE+00 

O-OOOOE+OO 

* 

0-7495E-23 

-0-2221E-20 

0.2076E-20 

0.3847E-03 

O-OOOOE+OO 

O-OOOOE+OO 

* 

O-OOOOE+OO 

O-OOCOE+OO 

O-OOOOE+OO 

O-COOOE+00 

0.2801E-03 

-0-9472E-10 

* 

O-OOOOE+GO 

O-OOOOE+OO 

O.QOOOE+00 

0.  GOOOE+00 

-0.9472E-10 

0-2S01E-C3 

*************************************************************************** 
******************  ********************************************************** 


* 

* 

El  = 

0-18E71 OD+06 

E2  = 

0-783349D+04 

E3 

0-783375D+04 

* 

* 

* 

NU23  = 

Q.42&386D+C0 

NU12  = 

0-297201  D+00 

NU13  = 

C-297186D+00 

IV 

* 

* 

* 

623  = 

0-259957D+C4 

Si  2  = 

0-357042D+04 

613  = 

C.357046D+04 

* 

* 

**************************************************************************** 


ROTATION  45  DEGREES  AROUND  XI 

******************************* 

**************************************************************************** 

*  HOMOGENIZED  COMPLIANCE  TENSOR  (BISECTING  DIRECTIONS)  * 


1  1 

2  2 

3  3 

2  3 

1  3 

1  2 

* 

0.529?£-05 

-0-1 575E-05 

-0-1575E-05 

0-8287E-1  0 

O-OOOOE+OO 

O-COOOE+OO 

* 

-0-1575E-05 

0-1 327E-03 

-0-5y69E-04 

-0-2108E-08 

O-COOOE+00 

O-OOOOE+OO 

* 

-0-1 575E-C5 

-0-5  969E-04 

0.1327E-03 

-0-2108E-08 

O-COOOE+OO 

O-OOOOE+OO 

* 

0-8287E-10 

-0-2103E-0i> 

-0-2108E-0b 

0-3647E-03 

O-OOOOE+OO 

O-OOOOE+OO 

* 

O-OOOOE+OO 

O-OOOOt+OO 

O-OOOOE+OO 

O-OOOOE+OO 

0-2301E-03 

-0-1869E-08 

* 

O-OOOOE+OO 

O-OOOOE+OO 

O-COOOE+00 

O-COOOE+00 

-0-1669E-08 

0-2801E-03 

* 

**************************************************************************** 

******************  ******************************  **************************** 
*  * 

•  El  =  0-13871 00+C6  E2  =  0. 7538400+04  E3  =  0-753840D+04  * 

*  * 

*  NU23  =  0-449934D+C0  NL12  =  0-2971 93D400  NU13  =  0-297193D+00  * 

•  * 

*  623  =  0-27421 1D+04  G1 2  =  0-3570440404  G13  =  0-3570440+04  * 

*  * 
************************************ **************************************** 
.GENERIC  NAME  : 

V3 


1 


I 


IF16  ? 


3 

lOPT  ? 


1 

2 

3 

A 

5 

6 

7 

8 


ISOTROPIC  CIRCULAR  FIBER) 

ORTHOTROPIC  CIRCULAR  FIBER) 

ISOTROPIC  KIDNEY  FIBER) 

ORTHOTROPIC  KIDNEY  FIBER) 

ISOTROPIC  STAGGERED  FIBER) 

ORTHOTROPIC  STAGGERED  FILER) 

HEXAGONAL  CELL  ISOTRCPIC  KIDNEY  FIBER) 
HEXAGONAL  CELL  ORTHOTROPIC  KIDNEY  FIBER) 


<IF  = 
CIF  = 


1  PROVIDE  RATIO  :  IMPREGNATED 

2  PROVIDE  RATIO  :  IMPREGNATED 

HDIST  :  MESH  DENSITY 


RESIN 

RESIN 

PARAMETER 


lOPT  =  ? 

1 

RATIO  =  ?  (  TXMIN  =  0.3545  ) 

.5 

E  ^  NU  =  ?  CYOUNG  MODULUS  POISSON  COEFF. 

84U00.  .22 

E  #  NU  =  ?  CYOUNG  MODULUS  POISSON  COEFF.  ISOTROPIC  RESIN) 
4000.  .34 


ISOTROPIC  FIBER) 


THE  FIBERS  ARE  PARALLEL  AT  XI 
*-***************«*  ************** 


NOTATIONS 
**  ******* 


cl ^  E  2#  c3 
GIJ 

NUCIJ) 

ER 

NUR 

NEL 

NOE 

TXR 


YOUNG  MODULI  FIBER  (COMPOSITE)  DIRECTIONS  X1^X2*>X 

FIBER  <  COMPOSITE  )  SHEAR  MODULI 

F13ER  <  COMPOSITE  )  POISSON  COEFFICIENT 

RESIN  YOUNG  MODULUS 

RESIN  POISSON  COEFFICIENT 

NUMBER  OF  ELEMENTS 

NUMBER  OF  NODES 

HcSiN  IMPREGNATED  RATIO  IN  VOLUME 


* 

ISOTROPIC  KIDNEY  (0.29)  FIBER 

•ISOTROPIC 

RESIN* 

MESH 

*. 

i*. 

* 

El  = 

84000.  .  612  =  34426.  .  NU12  - 

0.220 

« 

ER  = 

40C0.  * 

NEL  = 

310 

* 

£2  = 

84000-  .  G13  =  34426.  -  NU13  = 

0.220 

* 

NOE  = 

172 

* 

E3  = 

840C0.  .  G23  =  3^426.  .  NU23  = 

0.220 

* 

NUR  = 

0.340  * 

TXR  =  - 

5041 

*************************************************************************** 

*************************************************************************** 


*  HOMOGENIZED  ELASTIC  TENSOR 


* 

1  1 

2  2 

3  3 

2  3 

1  3 

1  2 

* 

* 

0.4673E+05 

0-5663E  +  04 

0-5663E+04 

0.3634E+03 

O.COOOE^OO 

O.OOOOE^OO 

* 

* 

0.5663E-t-04 

0.1 463E+05 

0.6551E^04 

0.9122E+03 

O.QOOOE-t-OO 

O.OOOOE<POO 

* 

* 

0.5663E^04 

0-6  55U+04 

0.l4o3E+05 

0.9122E+03 

O.OOOOEi-00 

O.OOOOE-POO 

* 

* 

0.3834C+03 

0.9122E+03 

0.9122E+03 

0-4503E+04 

O.OOOOE4>00 

O.OOOOE^OO 

* 

* 

O.OOOOE-t-OO 

O-OOCOE+00 

O-OOOOE+00 

O.COOOE+00 

0.4965E+04 

0.1195E+04 

* 

* 

O.OOOOE'>-O0 

O.OOOOE^OO 

O.COCOE-»-00 

O.OOOOE-t-OO 

0-1195E+04 

0.4965£'P04 

* 

REPRODUCED  AT  GOVERNMENT  EXPENSE  rO  ««*«««««« 


************* 

* 


HOMOGENIZED  COMPLIANCE  TENSOR 


*********************** 


* 

* 

* 

* 

* 

* 

* 

**************************************************** 


1  1 

0.2288E-04 

-0.6141E-05 

-0.6141E-05 

0.5396E-C6 

O-OOOOE+00 

0,0000c+00 


2  2  3  3  2  3 

-0.6141E-05  -0.6141E-05  0.5396E-06 

0,«734E-04  -0.3o20E-04  -0.9878E-05 
-0.3620E-04  0.S754E-04  -0.9«77E-05 

-0.9878E-05  -Q.9877E-05  0.2260E-03 

0,00C0E+00  O.COOOE+00  O.OOOOE+00 
O.OOOOE+OO  O.OOOOE+00  O.OOOOE+00  -0 


1  3 

.COOOE't'OO 
.OOOOE^OO 
.OOOOE+00 
.OOOOE+00 
.2133E-03 
.5146E-04 
****  ****  ** 


1  2 

O.OOOOE-i-00 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

-0.5146E-04 

0.213SE-03 

************* 


ORTHOTROPY  AXIS 

*************** 


HOMCGENI ZED 


COMPLIANCE  TEN 


************************ 

SOR  (BISECTING  DIRECTIO 


************** 
NS)  + 


********************** 

* 

*. 

* 

* 

* 

* 

* 

* 

* 

**************************************************************************** 


2  2 

3  3 

2  3 

1  3 

1  2 

* 

0.5371E-05 

-0.6411 E-05 

0.1930E-1  2 

O.OOOOE+OO 

O.OOOOE+OO 

* 

0.7230E-04 

-0.3083E-04 

0.36b9E-1 0 

O.OOOOE+OO 

O.OOOOE+OO 

* 

0.30S3E-04 

C.9205E-04 

-0.4395E-1 0 

O.OOOOE+OO 

O.OOOOE+OO 

* 

0.36fc9E-10 

-0.4595E-10 

0.2475E-03 

O.OOOOE+OO 

O.OOOOE+OO 

* 

O.OOOOE+OO 

O.QOOOE+00 

O.OOOOE+OO 

0.2652E-03 

-0.3174E-10 

* 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

-0.3174E-10 

0.1623E-03 

* 

1  1 

0.2283E-C4 
-0.587U-05 
-0.6411 £-05 
0.1930E-12 
O.OOOOE+OO 
O.OOOOE+OO 


**************************************************************************** 
*  * 

*  £1  =  0.436986D+05  E2  =  0.1383150+05  E3  =  C.108634D+05  * 

*  * 

+  NJ23  =  0.426484D+C0  N012  =  0.256556D+00  NU13  =  C.280135D+00  * 

*  * 

*  623  =  0.404050D+C4  612  =  0.o16080D+04  613  =  0.377019D+04  * 

*  * 
**************************************************************************** 

6ENERIC  NAME  : 

T4 

IFIB  ?  (  =  1  ISOTROPIC  CIRCULAR  FILER) 

<  =  2  ORTHOTROPIC  CIRCULAR  FIBER) 

(  =  3  ISOTROPIC  KIDNEY  FIBER) 

(  =  4  CRTHOTROPIC  KIDNEY  FIBER) 

(  =  5  ISOTROPIC  STAGGERED  FILER) 

(  =  6  CRTHOTROPIC  STAGGERED  FILER) 

(  =  7  HEXAGONAL  CELL  ISOTROPIC  KIDNEY  FIBER) 

C  =  b  HEXAGONAL  CELL  ORTHOTROPIC  KIDNEY  FIBER) 


4 

lOPT 


(IF  =  1  PROVIDE  RATIO 
(IF  =  2  PROVIDE  RATIO 

HOIST 

•> 


lOPT  = 

1 

RATIO  =  ?  (  TXMIN  =  0.3545  ) 

.5 


:  IMPREGNATED  RESIN 
:  IMPREGNATED  RESIN 
:  MESH  DENSITY  PARAMETER  ) 


asN34)A 


El  #  E2  #  E3  -  ?  (ORTHOTitOPIC  FIBERS  //  A  XI) 

380000.  14500.  14500. 

612  ^  613  ^  623  -  ?  (ORTHOTROPIC  FIBERS  H  A  XI) 

38300.  38000.  20000. 

NU12  #  NJ15  ^  N023  =  ?  (ORTHOTROPIC  FIBERS  //  A  XI) 

.22  .22  .25 

E  ^  NU  ?  (YOUNG  MODULUS  POISSON  COEFF.  ISOTROPIC  RESIN) 
3520.  .38 


THE  FIBERS  ARE  PARALLEL  AT  XI 
******************************** 


NOTATIONS  E1^E2^E3:  YOUNG  MODULI  FIBER  (COMPOSITE)  DIRECTIONS  X1#X2#X 

*********  GIJ  :  FIBER  (  COMPOSITE  )  SHEAR  MODULI 

NU(IJ)  :  FIBER  (  COMPOSITE  )  POISSON  COEFFICIENT 

ER  :  RESIN  YOUNG  MODULUS 

NUR  :  RESIN  POISSON  COEFFICIENT 

NEL  :  NUMBER  OF  ELEMENTS 

NOE  :  NUMBER  OF  NODES 

TXR  :  RESIN  IMPREGNATED  RATIO  IN  VOLUME 


*************************************************************************** 


* 

*  . 

KIDNEY  FIBER  (0-29) 

CRTHOTROPIC 

♦ISOTROPIC 

.* _ _ 

RESIN* 

MESH 

* 

El  = 

330000.  .  G12  =  38000. 

.  NU12  >  0.220 

« 

ER  = 

3520.  * 

NEL  * 

310 

* 

E2  = 

14500.  .  G13  =  33000. 

.  NU13  =  0.220 

* 

* 

NOE  = 

172 

* 

E3  = 

14500-  .  G23  =  2UG00. 

-  NU23  =  0-250 

* 

NUR  = 

0.330  * 

TXR  = 

.5041 

******************  **************************************  ******************* 


********************************************  ****  *************************** 


* 

HOMOGENIZED  ELASTIC  TENSOR 

* 

1  1 

2  2 

3  3 

2  3 

1  3 

1  2 

* 

0-1927E+06 

0.4129E-^04 

0-4129E+04 

0-1017E+02 

O.COOOE'i'OO 

0.0000E400 

* 

0.4129E+04 

0-9432E404 

0.4565E+04 

0.1699E^03 

O-CODOE-t-OO 

O-OOOOE+OO 

• 

0.4129E+04 

0-4  565E+04 

0.9432E+04 

0-1 899E+03 

0.Q000E400 

O.OOOOE'FOO 

* 

0.1017E+02 

0-1 899E+03 

0.1899E+03 

0.3536E+04 

O.COOOE-i'OO 

O.QOOOE400 

* 

0.0000E-«-00 

O.OOOOE+00 

O.COOOE+00 

O.OOOOE+OO 

0.4438E+04 

0.1145E+04 

* 

O.OOOOE-»-00 

O.OOOOE-^00 

O.OOOOE+00 

O.COOOE'i'OO 

0-1145E+04 

0.4438E-I-04 

************************************************************************** 


*************************************************************************** 

*  HOMOGENIZED  COMPLIANCE  TENSOR 


* 

1  1 

2  2 

3  3 

2  3 

1  3 

1  2 

* 

* 

0-5256E-05 

-0.1 553E-05 

-0-1553e-05 

0.1517E-06 

O.OOOOE+OO 

O.OOOOE+OO 

* 

* 

-0.1553E-05 

0.1 390£-03 

-0.6651E-04 

-0-3888E-05 

O.COOOE+OO 

O.OOOOE+OO 

* 

* 

-0.1553E-05 

-0.6651E-04 

0.1590E-03 

-0.3888E-05 

O.QOOOE+OO 

O.OOOOE+OO 

* 

* 

0.1517E-06 

-0-3d88E-05 

-0-3838E-05 

0-2832E-03 

O.COOOE+OO 

O.OOOOE+OO 

* 

* 

O.OOOOE-t-00 

O.OOOOE+OJ 

O.COOOE+OO 

O.OOOOE+OO 

0.241 4E-03 

-0.6225E-04 

* 

* 

O.OOOOE-t-00 

O.OOOOEt-OU 

O.CJOOL'f’OO 

O.OOOOE+OO 

-0-6225E-04 

0.2414E-03 

* 

ROOUCEO  AT  GOVERNMENT  EXPE^S^  «  *  *  «  «  «  «  «*«««**««*«  «« 


CRTHOTROPY  AXIS 


o 

o 

c 

n 


o 


*  HOMOSENIZED  COK 

*- .... .. ..........  ...... .. 


****** *«*«*«* **«***«*«  ************************** 

PLIANCE  TEKSOR  (BISECTING  DIRECTIONS) 


* 

1  1 

2  2 

3  5 

2  3 

1  3 

1  2 

* 

0.5256E-05 

-0.1 477E-05 

-0.1628E-05 

-0.1109E-1  5 

O.OOOOE-FOO 

O.OOOOE^OO 

* 

-0.1477E-05 

0.1 032E-03 

-0.345SE-04 

-0.1007E-12 

O.COOOE-t'OO 

O.OOOOE^OO 

* 

-0.1628E-05 

-0.3458E-04 

0.1109E-03 

0.1064E-12 

O.OOOOE-t-OO 

O.OOOOE'FOO 

* 

-0-1109E-15 

-0.1 007E-1 2 

G.1064E-12 

0.4110E-03 

O.OOOOE<(-O0 

O.OOOOE-FOO 

* 

O.OOOOE+00 

O.OOOOE'^OG 

O.OOOOE+00 

0.0000E4-00 

0.3036E-03 

0.4099E-12 

* 

o.oaooE<»-oo 

O.OGOOE+Oa 

O-COOOE+00 

0.O00OE4-00 

0.4099E-12 

0.1791E-03 

**************************************************************  ************** 


************************************************  **************************** 
*  * 


El  = 

0.190243D+G6 

E2  = 

0.969418D+04 

E3 

C.901470D+04 

* 

NU23  = 

0.3351 900+00 

NU12  = 

0.280952D+00 

NU13  = 

0.309803D+00 

* 

* 

623  = 

0.243537D+C4 

G12  s 

0.558297D+04 

G13  = 

0.329352D+04 

* 

*  * 


**************************************************************************** 
^GENERIC  NAME  : 

?T5 

IfIB  2  <  =  1  ISOTROPIC  CIRCULAR  FIBER) 

(  =  2  CRTriOTROPIC  CIRCULAR  FIBER) 

<  =  3  ISOTROPIC  KIDNEY  FIBER) 

(  =  4  ORTHOTRGPIC  KIDNEY  FIBER) 

(  =  5  ISOTROPIC  STAGGERED  FIBER) 

(  =  e  CkTHOTROPIC  STAGGERED  FIBER) 

(  =  7  HEXAGONAL  CELL  ISOTROPIC  KIDNEY  FIBER) 

^  (  =  a  HEXAGONAL  CELL  ORTHOTROPIC  KIDNEY  FIBER) 

lOPT  ?  (IF  =  1  PROVIDE  :  LENGTH  RATIO 

:  RESIN  INTEGRATED  RATIO 
(IF  =  2  PROVIDE  :  L  ^  RATIO  AND  HDIST  :  DENSITY 
^  NESH  PARAMETER  ) 

L  ? 

1.73 

RATIO  ?  (  >0.199  ) 

.5 

=  ?  (YOUNG  MODULUS  POISSON  COEFF.  ISOTROPIC  FIBER) 

B4000.  .22 

E^jr  NU  =  ?  (YOUNG  MODULUS  POISSON  COEFF.  ISOTROPIC  RESIN) 

4000.  .34 


THE  FIBERS  ARE  PARALLEL  AT  XI 

******************************** 


AT  GOVERNMENT  IXTEMSR  «  «  « 
*  ♦  t 


NOTATIONS 


EWE2#E3 

GIJ 

NU(IJ) 

ER 

NUR 

NEL 

NO£ 

TXR 


YOUNG  MODULI  FIBER  CCOMPOSITE)  DIRECTIONS  X1^X2#X3 

FIBER  (  COMPOSITE  )  SHEAR  MODULI 

FIBER  (  COMPOSITE  >  POISSON  COEFFICIENT 

RESIN  YOUNG  MODULUS 

RESIN  POISSON  COEFFICIENT 

NUMBER  OF  ELEMENTS 

NUMBER  OF  NODES 

RESIN  IMPREGNATED  RATIO  IN  VOLUME 


**************************************************************************** 


* 

*  • 

STAGGERED 

CIRCLE 

FIBER(1.75)  ISOTROPIC 

♦ISOTROPIC 

.  * _ 

RESIN* 

MESH 

« 

>♦ 

* 

El 

=  84000. 

.  G12 

=  34426.  .  NU12  a 

0.220 

* 

ER  = 

4000.  * 

NEL  = 

224 

* 

* 

E2 

=  84000. 

.  G13 

=  34426.  .  NU13  a 

0-220 

ik 

* 

NOE  = 

135 

* 

* 

E3 

=  84000. 

.  623 

a  34426.  .  NU23  = 

0.220 

* 

NUR  a 

0.340  ♦ 

TXR  = 

-5129 

* 

* 

*  • 

HOMOGENIZED  ELASTIC  TENSOR 

* 

•  ♦ 

* 

1  1 

2  2 

3  3 

2  3 

1  3 

1  2 

* 

* 

0-4566E+05 

0.5246E+04 

0.5259E+C4 

-0.8083E-1  3 

D.COOOE+00 

O.OOOOE+OO 

* 

* 

0.5246E+04 

0-1 352E+05 

0.5675E+04 

-0.1 6b7E-1 2 

O.COOOE+OO 

O.OOOOE+OO 

* 

* 

0.5259E+04 

0.5  675E+04 

0.1359E+05 

-0-2499E-1 2 

O.OOOOE+OO 

O.OOOOE+OO 

* 

* 

-0-80b3E-15 

-0.1687E-12 

-0.2499E-12 

0.3883E+04 

O.OOOOE+OO 

O.COOOE+OO 

* 

* 

O.OOOOE+OO 

O.OOOOE+OO 

C.COOOE+00 

O.OOOOE+OO 

0.4005E+04 

0.5464E-03 

* 

* 

O.OOOOE+OO 

O.OOCOE+00 

O.COOOE+OO 

O.COOOL+00 

0.5464E-03 

0.4016E+04 

* 

********************************************  ****  **********  ******  hit**  h*  It****  * 
****************************  ***************lt*  ***************  **************** 

*  HOMOGENIZED  COMPLIANCE  TENSOR  * 

*11  22  33  23  13  12  * 

*  0.2326E-04  -0.6359E-05  “0,e343E-05  -0.2005E-21  O.COOOE+OO  0,0Q00E*00  * 

*  -0.6359E-05  0.9140E-04  -0.35726-04  0.1539E-20  O.OOOOE+00  O.OOOOE+CO  * 

*  -0.€348E-05  -0.3572E-04  0,9098E-04  0.4172E-20  O.OOOOE+00  O.OOOOE+00  * 

*  -0.2005E-21  0.1539E-20  0.4172E-20  0.2576E-03  O.ODOOE*00  O.OOOOE+00  * 

*  O.OOOOc+00  O-OOOOE^OO  O.OJOOE+OO  O-OOOOE+00  0.2497E-03  -0-3397E-10  * 

*  O.OOOOE+OO  O-OOOOE+00  O.COOOE+OO  O.OOOOE+00  -0.3597E-10  0.2490E-03  * 

**************************************************************************** 


******************  ********************************************************1,1, 
*  * 

*  El  =  0.429907D+05  E2  =  0.1094140+05  E3  =  C.109914D+05  * 

*  * 

*  NU23  =  0.390776D+C0  NU12  =  0-2733660+00  NUI 3  =  C.272890D+00  * 

*  * 

*  G23  =  0.388259D+04  G1 2  =  0.4016180+04  613  =  0-4004900+04  * 

*  * 
**************************  **********************  **************************** 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


R3TAT1CN  45  DEGREES  AROUND  X1 

******************************* 


******************************************** *********«***«****A************ 

*  HOMOGENIZED  COMPLIANCE  TENSOR  (BISECTING  DIRECTIONS)  * 


1  1  2  2  3  3  2  3 

0-2326E-04  -0-6353E-05  -0-6353E-05  0.1152E-07 


0.0000E’<-00  O.OOOOE'i-00 


*  -0.6353E-C5  0-9213E-0A  -0.3665E-04  -0.2076E-06  O.OOOOE^OO  O.OOOOE+00 

*  -0.6353E-05  -0-3o65E-04  0-9213E-04  -0.2076E-06  O.OOOOE+00  O-OOOOE+00 

*  0.1152E-07  -0-2076E-06  -0.2076E-06  0.2538E-03  O.OOOOE+OO  O.OOOOE+00 

*  O-OOOOE+00  O-OOQOE+00  0.0000E-*-00  O.OOOOE-t-OD  0.2493E-03  0.3508E-06 

*  O-OOOOE+00  O.OOGOE+00  O.OOOOE+00  O.OOOOE+00  0.3508E-06  Q-2493E“03 

************************  *********************************«******************" 


**  ***«*«  **********  **************  ****************  **************************** 
*  * 

*  El  =  0.429907D+C5  E2  =  0.1085460+05  E3  =  0.108546D+05  * 

*  * 

*  NU23  =  0.397661D+C0  NUI 2  =  0.273138D+00  NUI 3  =  0.273138D+00  * 

*  * 

*  623  =  0.3)40000+04  612  =  0.4010530+04  Gl 3  =  0-401053D+04  * 

*  * 
******************  ********************************************************** 

GENERIC  NAME  : 

T6 

IFIB  ?  <  =  1  ISOTROPIC  CIRCULAR  FIbER) 

(  =  2  CRTHOTROPIC  CIRCULAR  FIBER) 

<  -  3  ISOTROPIC  KIDNEY  FIBER) 

<  *  4  ORTHOTROPIC  KIDNEY  FIBER) 

<  =  5  ISOTROPIC  STAGGERED  FIBER) 

(  =  6  CRTHOTROPIC  STAGGERED  FIBER) 

<  =  7  HEXAGONAL  CELL  ISOTROPIC  KIDNEY  FIBER) 

(  a  6  HEXmGONAL  CELL  ORTHOTROPIC  KIDNEY  FIBER) 

6 

lOPT  ?  (IF  =  1  PROVIDE  ;  LENGTH  RATIO 

:  RESIN  INTEGRATED  RATIO 
(IF  =  2  PROVIDE  i  L  ^  RATIO  AND  HDIST  :  DENSITY 
2  MESH  PARAMETER  ) 

L  ? 

1 

RATIO  ?  (  >0.342  ) 

.5 

HOIST  ?  (  <  0-15) 

.1 

El  ^  E2  ^  £3  =  ?  (ORTHOTRCPIC  FIBERS  //  A  Xl ) 

380000.  14500.  14500. 

612  ^  613  ^  62i  =  ?  CORTHOTROPIC  FIBERS  //  A  Xl ) 

36000.  38000.  20000. 

NU12  ^  NU13  ^  NU23  =  ?  (ORTHOTROPIC  FIBERS  //  A  XI) 

.22  .22  .25 

I  ^^NU  =  ?  (YOUNG  MODULUS  POISSON  COEFF.  ISOTROPIC  RESIN) 

3520.  .38 


THE  FIBERS  ARE  PARALLEL  AT  Xl 

******************************** 


REPRODUCED  AT  GOVERNMENT  EXPENSI 


NOTATIONS  E1^L2^E3:  YCUN3  MODULI  FIBER  {COMPOSITE)  DIRECTIONS  Xl^X2#X 

*********  GIJ  ;  FIBER  <  COMPOSITE  )  SHEAR  MODULI 

NU(IJ)  :  FIBER  (  COMPOSITE  )  POISSON  COEFFICIENT 

ER  :  RESIN  YOUNG  MODULUS 

NUR  :  RESIN  POISSON  COEFFICIENT 

NEL  :  NUMBER  OF  ELEMENTS 

NOE  :  NUMBER  OF  NODES 

TXR  ;  RESIN  IMPREGNATED  RATIO  IN  VOLUME 


*************************************************************************** 

*  STAGGERED  CIRCLE  FIBERCI.OO)  ORTHOTROPIC  ^ISOTROPIC  RESIN*  MAILLAGE 

*  E1  =  330000.  .  G12  =  38000.  .  NU12  =  0-220  *  ER  =  3520-  •  NEL  =  416 

*  E2  =  14500-  .  G13  =  38000.  .  NU13  =  0-220  *  *  NOE  =  233 

*  E3  =  14500-  .  G23  =  20COO.  .  NU23  =  0.250  *  NUR  =  0.380  *  TXR  =  .5083 

*************************************************************************** 


*  **  ***  «  **^********^******************************************* 


* 

1  1 

2  2 

3  3 

2  3 

1  3 

1  2 

ft 

* 

0-1911E-406 

0-4126E-404 

0.4126E-F04 

0-1870E-13 

O-COOOE+CO 

O.OOOOE+OO 

ft 

•k 

0.41  26E-404 

0.9131E-404 

0.4751E+04 

0.371 3E-1 3 

O.COOOE+OO 

O.OOOOE+OO 

ft 

* 

0.41 26E+04 

0.4751C-404 

0.9132E+04 

-0.6529E-1 4 

O.COUOE+00 

O.OOOOE+OO 

ft 

* 

0.1S70E-13 

0.3713E-13 

-0.t529E-14 

0.41256+04 

O.COOOE+OO 

O.OOOOE+OO 

ft 

-ft 

O.OOOOE+CO 

O.OOOOE+OO 

O.COOOE+00 

O.COOOE+00 

0.3562E+04 

0.1212E-03 

ft 

* 

O.OOOOE+00 

C.OOOOE+00 

C-COOOE-400 

O.COOOE+OO 

0.121 2E-03 

0.3562E+04 

ft 

*  HOMOGENIZED  COMPLIANCE  TENSOR  *  *******  ***  ****** 


1  1 
0.5300E-05 
-0.1 575E-05 
-0.1 575E-C5 
-0.1 235E-22 
O.OOOOE+00 
O.COOOE+OO 


*************** 


2  2 

-0.1 575E-05 
0.1 506E-03 
-0.7765E-04 
-0.1 472E-20 
O-OOOOE+00 
C.OOOOE+00 


3  3 

-0.1575E-05 

-0-77e5E-04 

0.1506E-C3 

0.V445E-21 

O.COCOE+00 

0.0300E+00 


2  3 

-0.1 235E-22 
-0-1472E-20 
0.9445E-21 
0.2424E-03 
O-OOOOE+00 
O.OOOOE+00 


1  3 

O.COOOE+00 

O-COOOE+00 

O.COOOE^OO 

O.COOOE+00 

0.2807E-03 

-0.9550E-11 


1  2 

O.OOOOE+00 

O-OOOOE+00 

O.OOOOE^OO 

O.OOOOE-^OO 

-0.9550E-11 

0.2S07E-03 


OOOE+OO  0.0300E+00  O.OOOOE+00  -0.9550E-11  0.2S07E-03 

********************************************************* 


ft 

ft 

El 

0.188674D+06 

ft 

NU23  = 

0.515522D+00 

ft 

ft 

G23  = 

0.41251 OD+04 

E2  =  0-663E94D+04 

NU12  =  0-297217D400 

G12  =  0.3561950404 


E3 

C-663952D+04 

ft 

ft 

A 

NU13  = 

C.297169D+00 

ft 

G13  = 

C.356230D+04 

W 

ft 

ft 

REPRODUCeD  AT  GOVERNMENT  EXPCNSM  ««**««« 


ROTATION  4S  DEGREES  AROUND  XI 

******************************* 


************************** 

*  HOMOGENIZED  COM 


************************************************ 

PLIANCE  TENSOR  (BISECTING  DIRECTIONS) 


1  1 

0.5300E-05 
-0.1 575E-05 
-0.1 575E-C5 
0.2549E-09 
O.OOOOE+00 
O.OOOOE-*-00 


2  2 

-0.1 575E-05 
0.9709E-04 
-0.2412E-04 
-0.6  483E-08 
O.OOOOE+00 
0.00Q0E-»-00 


3  3 

-0-1575E-05 

-0.2412E-04 

0.9709E-C4 

-0.6483E-03 

O.OOOOE+00 

O.QOOOE+00 


* 
* 
* 
* 
* 
* 


2  3 

0.254  9E-09 
-0.6483E-D8 
-0.6483E-0S 
0.4565E-03 
O.OOOOE+00 
O.OOOOE+00 


1  3 

O.OOOOE^OO 

O.COOOE+00 

O.OOOOE+00 

O.COOOE+00 

0.2S07E-03 

-0-1381E-07 


1  2 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.C000E4-00 

-0-1381E-07 

0.2807E-03 


************************************************************************** 


* 

* 

El 

0.1fc8674D+C6 

E2 

0.1029980*35 

E3 

0.102998D+05 

* 

NU23  = 

0.24S435D+00 

NU12  = 

0.2971930*00 

NU13  = 

0.297193D+00 

* 

* 

623  = 

0.219038D+C4 

G12  = 

0.3562130*04 

613  = 

0.356213D+04 

************ 
,  GENERIC  NAM 
^T7 


IFIB  ? 


=  ?  ( 


7 

RATIO 

.5 

£  ^  NU  =  ? 
84000.  .22 
E  NO  =  ? 
4000.  .34 


1 

2 

3 

4 

5 

6 
7 
£ 


************************************************************** 

c  « 

• 

ISOTROPIC  CIRCULAR  FI3ER) 

CRTHOTROPIC  CIRCULAR  FIBER) 

ISOTROPIC  KIDNEY  FI3ER) 

CRTHOTROPIC  KIDNEY  FIBER) 

ISOTROPIC  STAGGERED  FIBER) 

CRTHOTROPIC  STAGGERED  FIBER) 

HEXAGONAL  CELL  ISOTROPIC  KIDNEY  FIBER) 

HEXAGONAL  CELL  ORTHOTROPIC  KIDNEY  FIBER) 

TXKIN  =  0-3961  ) 

YOUNG  MODULUS  POISSON  COEFF.  ISOTROPIC  FIEER) 

YOUNG  MODULUS  POISSON  COEFF.  ISOTROPIC  RESIN) 


THE  FIBERS  ARE  PARALLEL  AT  XI 
******************************** 


NOTATIONS 

********* 


E1*E2^E3 

GIJ 

NU(IJ) 

ER 

NUR 

NEL 

NOE 

TXK 


YOUNG  MODULI  FIBER  (COMPOSITE)  DIRECTIONS  X1^X2*X3 

FIBER  (  COMPOSITE  )  SHEAR  MODULI 

FIBER  (  COMPOSITE  )  POISSON  COEFFICIENT 

RESIN  YOUNG  MODULUS 

RESIN  POISSON  COEFFICIENT 

NUMBER  OF  ELEMENTS 

NUMBER  OF  NODES 

RESIN  IMPREGNATED  RATIO  IN  VOLUME 


AT  GO\<pRNVE(JTi|E<P|NJjE,  ****«♦*♦* 

«  «  *  «  «  « 


*  « 


39 


m 

•v 

30 


*********************************************«*****************************<§ 


* 

*  < 

HEXAGONAL 

CELL 

FIBER  ISOTROPIC 

KIDNEY 

♦ISOTROPIC 

.  * _ 

RESIN* 

MESH 

* 

El  = 

84000. 

.  612 

=  34426.  .  NU12 

=  0.220 

ER  = 

4000.  * 

NEL  = 

498 

* 

E2  = 

84  000. 

.  G13 

—  34426.  .  NU13 

=  0.220 

* 

* 

NOE  = 

274 

* 

E3  = 

84000. 

.  623 

=  34426.  .  NU23 

=  0.220 

* 

NLR  = 

0.340  * 

TXR  = 

.5089 

* 

HOMOG 

lENIZEO  ELAS 

TIC  TENSOR 

* 

1  1 

2  2 

3  3 

2  3 

1  3 

1  2 

0.4630E+05 

0.5548E+04 

Q.554SE+04 

0.3217E-04 

O.COOOE+00 

O.OOOOE+OO 

* 

0.5546E+04 

0.1 449E+05 

0.6149E+04 

0.6630E-04 

O.COOOE+00 

O.OOOOE+OO 

* 

0.554aE+04 

0.6149E+04 

0.1449E+05 

0.9458E-04 

O.COOOE+OO 

O.OOOOE+OO 

* 

0.3217E-04 

0.6630E-04 

0.9453E-04 

0.4170E+04 

O.COOOE+00 

O.OOOOE+OO 

* 

O.OOOOE+00 

O.OOOOE+OO 

O.OOOOE+OO 

O.COOOE+00 

C.4507E+04 

0.1544E-02 

A 

O.OOOOE+00 

O.OOOOE+OQ 

Q.OOOOE+00 

O.COOOE+00 

0.1544E-02 

0.4507E+04 

* 

*  g 

HOMOGENIZED  COMPLIANCE  TcNSCR 

* 

m  * 

* 

1  1 

2  2 

3  3 

2  3 

1  3 

1  2 

* 

* 

0.2309E-04 

-0.62C7E-05 

-0.6207E-05 

0.6134E-1  3 

O.COOOE+OO 

O.OOOOE+OO 

* 

* 

-0.6207E-05 

0.8  5b5£-04 

-C.3406E-04 

-0.5446E-1 2 

C.OOOOE+00 

O.OOOOE+OO 

* 

* 

-0.6207E-05 

-0.3406E-04 

O.E5&5E-04 

-0.135&E-11 

D.COOOE+OO 

O.OOOOE+OO 

* 

* 

0.6134E-13 

-0.5446E-1 2 

-0.1353E-11 

0.2398E-03 

O.COOOE+OO 

O.COOOE+OO 

* 

* 

O-OOOOE+00 

O.OOCOE+OO 

G.OOOOE+CO 

O.OOOOL+OO 

0.2219E-03 

-0.7602E-10 

* 

* 

O.OOOOE+OO 

O.OOOOE+OO 

O.COOOE+OO 

O.COOOE+OO  - 

0.7602E-10 

0.2219E-03 

* 

************************** ************************************ ************** 
*  * 

*  £1  =  0.4331520+C5  E2  =  0.1164790+05  E3  =  0.116479D+05  * 

*  * 

*  NJ23  =  O.396722D  +  C0  NU12  =  0. 266842  O+OO  Kt;l3  =  0.2668420+00  + 

*  * 

*  623  =  0.4169720+04  612  =  0.4506930+04  613  =  0.4506930+04  * 

*  * 
**************************************************************************** 


*2I*IiS**i5*S^S!*l****S2*2***’^2 


*************************************************************************** 

HOMOGENIItO  COMPLIANCE  TENSOR  (BISECTING  DIRECTIONS)  * 

________  ____  __  _  A 


A^GOVfRN^EN^EXP^NS^  ****♦♦**«♦ 


*  0.2309E-04  -0.6207E-05  -0.6207E-05  -0.4191E-13  O-OOOOE^OO  O-OOOOE+00 

*  -0.6207E-05  0.8565E-04  -0.3406E-04  -0.2220E-12  O-OOOOE^OO  0-0000E*00 

*  -0.6207E-05  -0.3405E-04  0.E585E-04  0-5913E-12  O-OOOOE+00  O-OOOOE+OO 

*  -0.4191E-13  -0.2220E-12  0.5913E-12  C.2398E-03  0-C000E>»-00  O-OOOOE+00 

*  O-OOOOE+00  O.OOOOE+00  O-COOOE+00  O.OOOOE+00  0.2219E-03  -0.3729E-10 

*  O-OOOOE+CO  O.OOOOE+OO  O.OOOOE+00  O.OOOOE+00  -0.3729E-10  0-2219E-03 

************************ **************************** «**«***«  ************* 


**  ******************************************  ****  *************************** 
* 


El 

0.4331 52D+05 

£2  = 

0.116479D+05 

E3 

C.116479D405 

NU23  = 

0-39672 2D+C0 

NU12  = 

0.2688420400 

NU13  = 

0.2688420400 

623  = 

0-416972D't-04 

612  = 

0.450693D404 

G13  = 

0.4S0692D404 

* 

******************************  ********************************************* 
,  CENERIC  NAME  : 

IFIB  ?  (  =  1  ISOTROPIC  CIRCULAR  FIBER) 

(  =  2  ORTriOTROPIC  CIRCULAR  FIBER) 

<  =  3  ISOTROPIC  KIDNEY  FIBER) 

(  =  4  CRTHOTROPIC  KIDNEY  FIBER) 

<  =  5  ISOTROPIC  STAGGERED  FIBER) 

(  =  6  CRTHOTROPIC  STAGGERED  FIBER) 

(  =  7  HEXAGONAL  CELL  ISOTROPIC  KIDNEY  FIBER) 

C  s  8  HEXAGONAL  CELL  ORTHOTROPIC  KIDNEY  FIBER) 

6 

RATIO  =  ?  (  TXHIN  =  0.3961  ) 

.5 

El  #  E2  /  E3  =  ?  (ORTHOTROPIC  FIBERS  //  A  XI) 

360000.  1a500.  1450C. 

612  #  613  G25  =  ?  (ORTHOTROPIC  FIBERS  //  A  X1> 

36000.  38C00.  20000. 

NU12  #  NU13  *  NU23  -  ?  (ORTHOTROPIC  FIBERS  t!  A  XI) 

.22  .22  .25 

E  #  NU  =  ?  (Y0UN3  MODULUS  POISSON  COEFF-  ISOTROPIC  RESIN) 

3520.  .38 


12 f**J *121***21*** *2***; ***1**^1 


NOTATIONS 

**  ******* 


E1^E2^E3:  YOUNG  MODULI  FIBER  (COMPOSITE)  DIRECTIONS  X1*X2*X3 

GIJ  :  FIBER  (  COMPOSITE  )  SHEAR  MCDUH 

NU(IJ)  :  FIBER  (  COMPOSITE  )  POISSON  COEFFICIENT 

ER  :  RESIN  YOUNG  MODULUS 

NUR  :  RESIN  POISSON  COEFFICIENT 

NEL  ;  NUMBER  OF  ELEMENTS 

NOE  :  NUMBER  OF  NODES 

TXR  :  RESIN  IMPREGNATED  RATIO  IN  VOLUME 


******************  ***«***«*-^*«  ************************  ********************** 

•HEXAGONAL  CELL  FI6ER  ORTHOTROPIC  KIDNEY  *ISOTROPIC  RESIN*  MESH  * 


REPRqpuqjtV^qPYERNMiNX  «  « 


»  »  »  «■  «  »  »  »  »»»»»»»*'»»»  »  »  »  »  *  » 


*  tl  =  3tiC000.  .  (>12  =  3o000. 

*  E2  >  14500.  .  G13  «  33000. 

*  E3  >  14500.  .  G23  «  20C00. 


NL12  = 

0.220 

« 

ER  = 

3520. 

NEL  = 

498  + 

NU13  = 

0.220 

* 

« 

NOE  * 

274  * 

NU23  = 

0.250 

* 

NUR  = 

0.380 

« 

TXR  » 

.5089  * 

X 

* 


NONCGENIZEO  ELASTIC  TEKSOR 


1  1 

0.1909k+06 

0.4127E+04 

0.4l2ttE+04 

-0.2039E-05 

O.OOOOE+00 

O.OOOOE+00 


2  2 

0.4127E+04 
0.9  448E  +  04 
C.4488E-f04 
0.7904E-05 
O.OOOOE+00 
O.OOOOE+OO 


3  3  2  3 

0.4123E+04  -0.2039E-05 
0.448SE+04  0.7904E-05 

0-9477C+04  0.2161E-04 

C.2161E-04  0.3341E+04 

O.OOOOE+00  O.OOOOE-^OO 
O.COOOE+00  O.OOOOE+OO 


1  3 

O.OOOOE+OO 

O.COOOE+00 

O.OOOOE+OO 

O.COOOE+00 

0.3980E+04 

0.2958E-02 


1  2 

O.COOOE+00 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.2955E-02 

0.3980E+04 


o 

*  n 

.*0 

* 

* 

* 

* 

* 

* 

* 


*  **  *  ^ 

.* 
* 
* 


HOMOGENIZED  COMPLIANCE  TENSOR 


1  3 
O.OOOOE+OO 
O.OOOOt+00 
O.OOOOE+OO 
O.OOOOE+OO 


1  2 
O.OOOOE+OO 
O.OOOOE+OO 
O.OOOOE+OO 
O.OOOOE+OO 


* 

* 

* 

* 


1  1  2  2  3  3  2  3 

0.5307E-05  -0.1574E-05  -0.1566E-05  0.1710E-13 

-0.1574E-05  0.1370E-03  -0.6420E-04  0.9021E-13 

-0.1566E-05  “0.6420E-04  0.1366E-03  -0.7328E“12 

0.1710E-13  0.9021E-13  -0.732oE-12  0.2993E-03  .  ,  „ 

O.OOOOE+OO  O.OOOOE+OO  O.COOOE+00  O.COOOE+DO  0.2513E-03  -0.1867E-09 

O.OOOOE+OO  O.OOOOE+OO  O.COCOE+00  O.COOOt+OO  -0.1867E-09  0.2513E-03 

*************************************************************************** 

*************************************************************************** 

* 

El  =  0.138440D+C6  E2  =  0.729813D+04  E3  =  C.732080D+04  * 

* 
* 
* 
* 
* 


NJ23  =  0.46851 OD+CO 

623  =  0.334081D+C4 


NU12  =  0.2966740+00 

G12  =  0. 3979740+04 


NU13  =  C. 2950990+00 
G13  =  0.3979740+04 


**********  ****************************************************************** 


ROTATION  45  DEGREES  AROUND  X1 

******************************* 

**************************************************************************** 


*  HCMCGENIZEO  COMPLIANCE  TENSOR  (BISECTING  DIRECTIONS) 


*. 

* 

1  1 

2  2 

3  3  '■ 

2  3 

1  3 

1  2 

* 

* 

0.53076-05 

-0.1 570E-05 

-0.1570E-05 

0.8355E-08 

O.COOOE+00 

O.OOOOE+OO 

* 

* 

-0.1 570E-05 

0.1  111 E-03 

-0.38536-04 

-0.2122E-06 

O.OOOOE+OO 

O.OOOOE+OO 

* 

* 

-0.1 570E-05 

-0.3853E-04 

0.1111E-03 

-0.2122E-06 

O.OOOOE+OO 

O.COOOE+00 

* 

0.8355E-08 

-0.2122E-06 

-0.2122E-06 

0.4020E-03 

O.COOOE+00 

O.COOOE+OO 

* 

* 

O.COOOE+00 

O.OOOOE+OO 

O.OOOUE+00 

O.OOOOE+OO 

0.251 3E-03 

-0.1117E-10 

* 

* 

O.CUOOE+00 

O.OOOOE+OO 

O.COOOE+Od 

O.OOCOt+00 

-0.1117E-10 

0.2513E-03 

* 

* 

* 

« 

* 

* 

* 

* 

* 


El 

0.1884400+06 

E2  = 

0.6997750+04 

NU23  * 

0.3466430+00 

NU12  = 

0.2958860+00 

G23  = 

0.2487500+04 

G12  = 

0.3979740+04 

**************************** 

▲ 

E3 

0.8997750+04 

* 

NU13  = 

C. 2958860+00 

* 

* 

G13  = 

0.3979730+04 

* 

* 

AT  GOVERNMENT  EXPENSE 


I 


Listing  2:  Evaluation  of  Microstress  Fields 


C  0  If  P  0  15-JUL-1087  at  12:21:55 


HOMOGENIZATION  LAMINATED  COMPOSITE 


******************************************************************************* 


*  STAGGERED  FIBERS  PARALLEL  TO  THE  AXIS  X3  * 

* 

*  N  N  * 

*  (  X  )  (  Y  )  A  =  0.282094790-03  * 

*  FIBER  SECTION  :  ( - )  +  ( - )  =  1  ;  WITH  :  B  =  0.28209479D-03  * 

*  (  A  )  (  B  )  N  =  0.20000000D+01  * 

*  * 

*  SPACING  OF  FIBERS  :  IN  X  =  O.lOOOD-02  IN  Y  =  O.lOOOD-02  ♦ 

♦ 

«  NUMBER  OF  ELEMENTS  ^  104  NUMBER  OF  NODES  =  69  RESIN  RATIO  =  .5130  * 

*  ♦ 

♦ . * 

*  ISOTROPIC  RESIN  CHARACTERISTICS  * 

*  * 

♦  YOUNG’S  MODULUS  =  3520.0  POISSON’S  COEFFICIENT  =  0.38000  ♦ 

♦  * 

* . . . * 

*  ORTHOTROPIC  FIBER  CHARACTERISTICS  * 

♦  * 

*  El  =  0.380000D+06  E2  =  0.145000D+05  E3  »  0.145000D+05  * 

*  G23  =  0.200000D+05  G13  =  0.380000D+05  G12  ==  0.380000D+05  * 

*  NU23  =  0.250000D+00  NU13  =  0.220000D+00  NU12  =  0.220000D+00  * 

*  * 
*^,*^,i,4,*^i,********************************************************************** 

4****************************************************************************** 

*  HOMOGENIZED  ELASTIC  TENSOR  * 

*  . . * 

*  11  22  33  23  13  12  * 

*  0.1336E+05  0.5647E+04  0.4021E+04  O.OOOOE+00  O.OOOOE+00  0.3978E-12  * 

«  0.0492E+04  0.399gE404  O.OOOOE-i-00  O.OOOOE-^00  -0.5453E-13 

*  0.1094E+05  O.OOOOE+00  O.OOOOE+00  0.6360E-14  * 

*  0.3472E+04  0.3325E-03  O.OOOOE+00  * 

*  0.3767E+04  O.OOOOE+00  * 

*  0.5600E+04  * 

:t****:Hi;H,****t*lH,****t********************************************************* 

*t*0t^t*^1,**4,t01tm*************************************************************** 

*  HOMOGENIZED  COMPLIANCE  TENSOR  * 

*  . . . ....* 

*  11  22  33  23  13  12  * 

*  0.1030E-03  -0.5350E-04  -0.1828E-04  O.OOOOE-i^OO  O.OOOOE+00  -0.7823E-20  * 

*  0.1524E-03  -0.3604E-04  O.OOOOE-^00  O.OOOOE-fOO  0.5326E-20  o 

>»  0.1113E-03  O.OOOOE+00  O.OOOOE+OO  0.840gE-21  * 

*  0.2880E-03  -0.2542E-10  O.OOOOE-^00  * 

*  0.2655E-03  O.OOOOE^OO  * 

*  0.1788E-03  * 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


0tl***************************************************************************** 


i^^ti^^^mt^^^^tL^^^iti*************************************************************** 

N  EQUIVALENT  MODULI  • 


El  0.070682D+04 
G23  =  0.347231D+04 
NU23  ==  0.2364480-^00 


E2  =  0.6560660-1^04 
G13  =  0.3767000+04 
NU13  =  0.1774220+00 


* 

E3  =  0.8981110+04  * 
G12  =  0.6599010+04  * 
NU12  =  0. 6201840+00  ♦ 


^t**t^t^t************************************************************************ 
!  ! 

!  MESH  ! 

I  ! 

4i«4i:»4i4i4i*4<4>*4i**4i4i4>***4i*4i4>*  +  +  +  +  +  +  +  +  +  +  +  +  *  +  *  +  ***  +  +  **  +  **  +  +  +  +  «  +  + +  +  +  +  +  +  +  «<  +  +  +  +  +  +  +  +  +  +  +  ** 


number  of  elements  :  104  NUMBER  OF  NODES  :  28 


SO  =  1  — >  FIBER 
SO  =  2  — >  RESIN 


RTiE 

1 

SO 

1 

NODES 

1 

NODES  COORDINATES 

1 

1 

1 

! 

2 

1 

3 

0 

! 

X 

:-5.000E-04  -5.000E-04  -3.590E-04 

O.OOOE+00 

! 

Y 

;-3.S90E-04  -5.000B-04  ~5.000E-04 

O.OOOE+00 

2 

! 

1 

j 

2 

3 

10 

0 

1 

X 

:-S.OOOE-04  -3.590E-04  -3.005B-04 

O.OOOE+00 

i 

y 

;-3.590E-04  -5.000E-04  -3.005E-04 

O.OOOE+00 

3 

j 

1 

1 

2 

9 

8 

0 

! 

X 

:-S.000E-04  -3.839E-04  -5.000E-04 

O.OOOB+00 

i 

Y 

:-3.590E-04  -2.429E-04  -2.179E-04 

O.OOOE+00 

4 

! 

1 

! 

2 

10 

9 

0 

j 

X 

:-5.000E-04  -3.005E-04  -3.839E-04 

O.OOOE+OO 

! 

Y 

:-3.590E-04  -3.005E-04  -2.429E-04 

O.OOOE+OO 

5 

! 

1 

! 

10 

3 
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1  8.30 

! 

1 

37 

!  -74.3 

i 

104. 

!  -26.9 

!  O.OOOE+00 

O.OOOE+00 

1  -4.56 

! 

» 

38 

!  53.1 

! 

96.8 

!  44.3 

!  O.OOOE+00 

O.OOOE+00 

1  -2.73 

j 

1 

39 

!  24.2 

! 

106. 

!  36.9 

!  O.OOOE+00 

O.OOOE+00 

1  0.737 

! 

! 

40 

!  10.2 

i 

102. 

!  29.8 

!  O.OOOE+00 

O.OOOE+00 

1  -1.52 

j 

1 

41 

!  -27.9 

1 

66.5 

!  1.98 

!  O.OOOE+00 

O.OOOE+00 

1  -11.7 

i 

i 

42 

•  23.9 

i 

119. 

!  41.8 

?  O.OOOE+00 

O.OOOE+00 

1  0.852 

1 

! 

43 

!  -12.0 

! 

89.6 

!  -30.0 

!  O.OOOE+00 

O.OOOE+00 

1  -26.0 

1 

1 

44 

!  -74.3 

1 

104. 

!  -26.9 

!  O.OOOE+00 

O.OOOE+00 

1  4.56 

1 

1 

45 

!  -15.5 

1 

117. 

!  26.0 

!  O.OOOE+00 

O.OOOE+00 

!  4.683E-15 

1 

! 

46 

!  -2.12 

! 

81.3 

!  17.4 

!  O.OOOE+00 

O.OOOE+00 

1  -2.71 

1 

1 

47 

!  -8.75 

i 

26.2 

!  -6.06 

!  O.OOOE+00 

O.OOOE+00 

1  1.29 

1 

! 

48 

!  -11.2 

! 

126. 

!  -20.9 

!  O.OOOE+00 

O.OOOE+00 

1  0.994 

1 

j 

49 

!  37.5 

i 

112. 

!  -24.0 

!  O.OOOE+00 

O.OOOE+00 

1  1.583E-14 

1 

! 

50 

!  -9.76 

! 

117. 

!  -23.2 

!  O.OOOE+00 

!  O.OOOE+00 

1  -8.30 

1 

i 

51 

!  37.5 

! 

112. 

!  -24.0 

!  O.OOOE+00 

!  O.OOOE+00 

1-1.446E-14 

1 

I 

52 

!  39.0 

! 

116. 

!  -23.0 

!  O.OOOE+00 

!  O.OOOE+00 

1  -19.8 

1 

1 

53 

!  -9.76 

! 

117. 

1  -23.2 

!  O.OOOE+00 

!  O.OOOE+00 

1  -8.30 

1 

! 

54 

!  37.5 

! 

112. 

I  -24.0 

!  O.OOOE+00 

I  O.OOOE+00 

1-1.289E-14 

1 

! 

55 

!  -15.6 

! 

117. 

1  26.0 

!  O.OOOE+00 

!  O.OOOE+00 

1-1.878E-15 

1 

! 

56 

f  -2.12 

I 

81.3 

f  17.4 

!  O.OOOE+00 

!  O.OOOE+00 

1  -2.71 

1 

! 

57 

!  -11.2 

! 

126. 

!  -20.9 

1  O.OOOE+00 

1  O.OOOE+00 

1  0.994 

1 

1 

58 

!  -8.75 

1 

26.2 

1  -6.06 

!  O.OOOE+00 

1  O.OOOE+00 

1  1.29 

1 

! 

50 

!  30.0 

!  116. 

!  -23.0 

!  O.OOOE+00 

! 

O.OOOE+00 

!  -19.8 

1 

60 

!  37.5 

I  112. 

!  -24.0 

1  O.OOOE+00 

! 

O.OOOE+00 

!  1.739E-14 

! 

61 

1  53.1 

!  06.8 

1  44.3 

!  O.OOOE+00 

! 

O.OOOE+00 

!  -2.73 

! 

62 

1  10.2 

1  102. 

1  29.8 

!  O.OOOE+00 

! 

O.OOOE+00 

!  -1.52 

! 

63 

!  -27.9 

i  66.5 

!  1.08 

!  O.OOOE+00 

! 

O.OOOE+00 

!  -11.7 

f 

64 

!  24.2 

!  106. 

!  36.0 

1  O.OOOE+00 

I 

O.OOOE+00 

!  0.737 

! 

65 

!  23.0 

1  110. 

!  41.8 

!  O.OOOE+00 

! 

O.OOOE+00 

1  0.852 

! 

66 

!  -12.0 

!  89.6 

!  -30.0 

f  O.OOOE+00 

1 

O.OOOE+00 

1  -26.0 

! 

67 

!  -9.76 

!  117. 

!  -23.2 

!  O.OOOE+00 

1 

O.OOOE+00 

i  8.30 

1 

68 

!  -74.3 

!  104. 

!  -26.0 

f  O.OOOE+00 

! 

O.OOOE+00 

!  -4.56 

i 

69 

!  -12.0 

!  89.6 

!  -30.0 

!  O.OOOE+00 

! 

O.OOOE+00 

!  26.0 

1 

70 

!  -74.3 

!  104. 

!  -26.0 

!  O.OOOE+00 

1 

O.OOOE+00 

!  4.56 

! 

71 

!  39.0 

!  116. 

!  -23.0 

!  O.OOOE+00 

! 

O.OOOE+00 

!  10.8 

i 

72 

!  -11.2 

!  126. 

!  -20.9 

!  O.OOOE+00 

! 

O.OOOE+00 

! -0.994 

! 

73 

!  -15.5 

!  117. 

!  26.0 

!  O.OOOE+00 

I 

O.OOOE+00 

!  4.93gE-15 

1 

74 

!  -2.12 

!  81.3 

!  17.4 

!  O.OOOE+00 

j 

O.OOOE+00 

!  2.71 

1 

75 

!  -8.75 

!  26.2 

!  -6.06 

!  O.OOOE+00 

1 

O.OOOE+00 

!  -1.29 

! 

76 

!  -39.7 

!  106. 

!  -26.2 

!  O.OOOE+00 

I 

O.OOOE+00 

!  -5.81 

1 

77 

!  -17.3 

!  90.8 

!  -27.5 

!  O.OOOE+00 

J 

O.OOOE+00 

!  -18.8 

! 

78 

!  19.2 

!  118. 

!  -22.6 

!  O.OOOE+00 

1 

O.OOOE+00 

!  -3.73 

! 

79 

!  -15.5 

!  66.1 

!  6.52 

!  O.OOOE+00 

! 

O.OOOE+00 

!  -10.1 

i 

80 

!  -11.2 

!  117. 

!  -23.0 

!  O.OOOE+00 

! 

O.OOOE+00 

!  -5.72 

! 

81 

!  3.42 

!  92.0 

!  23.9 

!  O.OOOE+00 

1 

O.OOOE+00 

!  -4.59 

! 

82 

!  8.38 

!  95.3 

!  26.7 

!  O.OOOE+00 

1 

O.OOOE+00 

!  1.62 

! 

1 

83 

!  58.0 

!  105. 

!  49.4 

!  O.OOOE+00 

! 

O.OOOE+00 

!  -2.78 

! 

84 

!  58.0 

!  105. 

!  49.4 

!  O.OOOE+00 

I 

O.OOOE+00 

!  2.78 

! 

! 

85 

!  -27.9 

!  66.5 

!  1.98 

!  O.OOOE+00 

j 

O.OOOE+00 

!  11.7 

! 

i 

86 

!  23.0 

!  119. 

!  41.8 

!  O.OOOE+00 

j 

O.OOOE+00 

1-0.852 

j 

87 

J  3.42 

1  92.9 

!  23.9 

!  O.OOOE+00 

j 

O.OOOE+00 

1  4.59 

! 

! 

88 

!  8.38 

!  95.3 

!  26,7 

!  O.OOOE+00 

! 

O.OOOE+00 

!  -1.62 

! 

1 

89 

!  -15.5 

!  66.1 

!  6.52 

}  O.OOOE+00 

} 

O.OOOE+00 

1  10.1 

I 

! 

90 

!  24.2 

!  106. 

!  36.9 

!  O.OOOB+00 

1 

O.OOOE+00 

1-0.737 

! 

! 

91 

!  10.2 

!  102. 

!  29.8 

!  O.OOOE+00 

! 

O.OOOE+00 

1  1.52 

1 

92 

!  53.1 

!  06.8 

!  44.3 

!  O.OOOE+00 

1 

O.OOOE+00 

1  2.73 

1 

! 

93 

!  16.2 

!  114. 

!  -23.6 

!  O.OOOE+00 

j 

O.OOOE+00 

1-4.016E-14 

J 

1 

94 

!  -9.26 

1  24.5 

!  -6.90 

!  O.OOOE+00 

! 

O.OOOE+00 

1  1.24 

J 

1 

05 

!  -11.3 

!  121. 

!  -22.2 

!  O.OOOE+00 

t 

O.OOOE+00 

1  0.590 

1 

96 

!  -20.1 

!  110. 

!  21.6 

f  O.OOOE+00 

j 

O.OOOE+00 

1  6.640E-15 

! 

1 

97 

!  -20,1 

!  110. 

!  21.6 

!  O.OOOE+00 

1 

O.OOOE+00 

1  5.927E-15 

! 

98 

!  -17.3 

!  09.8 

!  -27.5 

!  O.OOOE+00 

! 

O.OOOE+00 

1  18.8 

! 

1 

09 

!  -11.2 

!  117. 

!  -23.0 

!  O.OOOE+00 

1 

O.OOOE+00 

1  5.72 

1 

! 

100 

!  19.2 

!  118. 

!  -22.6 

!  O.OOOE+00 

\ 

O.OOOE+00 

1  3.73 

! 

1 

101 

!  -9.26 

!  24.5 

!  -6.90 

!  O.OOOE+00 

! 

O.OOOE+00 

1  -1.24 

1 

102 

!  -11.3 

!  121 . 

!  -22.2 

!  O.OOOE+00 

j 

O.OOOE+00 

1-0.590 

! 

103 

!  -39.7 

!  106. 

!  -26.2 

!  O.OOOE+00 

! 

O.OOOE+00 

1  5.81 

! 

104 

!  16.2 

!  114. 

!  -23.6 

!  O.OOOE+00 

! 

O.OOOE+00 

1-2.327E-14 

• 

1 

! 

1 

STRESS 

FORCE 

1 

1 

1. 

NODE  !  X 

! 

Y  ! 

FX  ! 

FY  ! 

FZ 

1  !  -3.8394E-04  !  -2. 

428gE-04  1 

1.015  ! 

-44.57  ! 

O.OOOOE+00 

2  !  -3.0063E-04  !  -3. 

0053E-04  ! 

16.65  1 

-52.27  ! 

O.OOOOE+00 

3  !  -2.4289E-04  !  -3. 

8304E-O4  1  - 

B.023  ! 

-24.66  1 

O.OOOOEh^OO 

4  !  2.4289E-04  1  -3.83g4£-04  1 

8.023  ! 

-24.66  ! 

O.OOOOE+00 

5  I  -2.6711E-04  !  -1. 

1606E-04  1 

14.73  ! 

33.01  1 

C  OOOOE+00 

6  !  -1.9947E-04  1  -1. 

9947E-04  1  - 

16.87  1 

52.66  ! 

O.OOOOE+00 

7  !  -1.1606E-04  1  -2. 

6711E-04  !  - 

10.48  ! 

35.70  ! 

O.OOOOE-eOO 

8  1  3.8304E-04  !  -2. 

4280E-04  1  - 

1.915  I 

-44.57  I 

O.OOOOE+00 
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9 

! 

3.0053E-04 

1 

-3.0053E-04 

1 

-16.66 

!  -62.27 

!  O.OOOOE+00 

I 

10 

! 

1.1606E-04 

! 

-2.6711E-04 

! 

10.46 

!  35.70 

1  O.OOOOE+00 

1 

1  • 

11 

1 

O.COOOE+00 

! 

-2.8209E-04 

1 

-4.7635E-06 

!  63.08 

t  O.OOOOEh^OO 

I 

12 

! 

-2.8200E-04 

! 

O.OOOOE+00 

! 

-6.002 

i  1.8324E-07 

1  O.OOOOE-t-00 

! 

13 

1 

2.6711E-04 

1 

-1.1606E-04 

( 

-14.73 

t  33.91 

1  O.OOOOE+00 

1 

14 

! 

1.0047E-04 

! 

-1.9947E-04 

! 

16.87 

1  62.66 

I  O.OOOOE+00 

! 

16 

1 

-2.6711E-04 

1 

1.1606E-04 

! 

14.73 

!  -33.91 

1  O.OOOOE-pOO 

I 

16 

! 

-1.0g47E-04 

! 

1.0047E-04 

! 

-16.87 

!  -62.66 

I  O.OOOOE-pOO 

! 

17 

! 

2.820gE-04 

! 

O.OOOOE+00 

! 

6.002 

!  -1.8324E-07 

1  O.OOOOE+OO 

! 

18 

! 

-3.8394E-04 

1 

2.4289E-04 

! 

1.015 

]  44.67 

!  O.OOOOE+00 

1 

10 

! 

-3.00S3E-04 

i 

3.0053E-04 

! 

16.66 

1  62.27 

!  O.OOOOE+OO 

! 

20 

i 

-1.1606E-04 

! 

2.5711E-04 

1 

-10.46 

!  -35.70 

!  O.OOOOE+00 

! 

21 

j 

O.OOOOE+00 

1 

2.820gE-04 

! 

4.7635B-06 

!  -63.08 

!  O.OOOOE+00 

j 

! 

22 

! 

1.9947E-04 

i 

1.9947E-04 

i 

16.87 

!  -62.66 

!  O.OOOOE+OO 

! 

1 

23 

! 

1.1606E-04 

1 

2.6711E-04 

! 

10.46 

!  -35.70 

!  O.OOOOE-pOO 

j 

24 

1 

2.5711E-04 

! 

1.1606E-04 

! 

-14.73 

!  -33.91 

!  O.OOOOE+00 

J 

25 

! 

-2.428gE-04 

I 

3.83g4E-04 

I 

-8.023 

!  24.66 

!  O.OOOOE+00 

j 

26 

! 

3.0053E-04 

! 

3.0053E-04 

1 

-16.65 

!  52.27 

!  O.OOOOE+00 

j 

27 

1 

2.428gE-04 

! 

3.83g4E-04 

1 

8.023 

!  24.66 

!  O.OOOOE+00 

i 

28 

i 

3.8394E-04 

! 

2.428gE-04 

i 

-1.916 

!  44.67 

!  O.OOOOE^^OO 

! 

^lUttif^l^Li^^^^l^^^if^L^tL************************************************************** 

!  ! 


MACROSCOPIC  STRESS  FIELD 


! 


S(l,l)  =  O.OOOOD+00  S(l,2)  =  O.OOOOD+00  S(l,3)  =  O.lOOOD+03  ! 

S(2,2)  =  O.OOOOD+00  S(2,3)  =  O.OOOOD+OO  ! 

S(3,3)  =  O.OOOOD+00  ! 


! 


MACROSCOPIC  STRAIN  TENSOR 


! 

I 


E(l,l)  =  -0.2436D-02  E(l,2)  =  -0.3604D-02  E(l,3)  =  0.1113D-01  ! 

E(2,2)  =  O.OOOOD+00  E(2,3)  =  O.OOOOD+00  ! 

E(3,3)  =  O.OOOOD+00  ! 

_ ! 

i 

MICROSCOPIC  STRESS  FIELD  ! 

! 


! 

ELE 

Sll 

!  S22 

!  S33 

i  S23 

! 

S13 

! 

S12 

! 

1 

1 

1 

-6.35 

!  -7.92 

!  150. 

1  O.OOOE+00 

! 

O.OOOE-pOO 

1 

6.210E-16  i 

! 

2 

-19.8 

!  -8.08 

!  159. 

!  O.OOOE+00 

1 

O.OOOE+00 

1 

-1.35 

! 

3 

-26.4 

!  -3.32 

!  160. 

!  O.OOOE+00 

i 

O.OOOE+OO 

! 

-10.4 

! 

4 

-21.8 

!  -7.27 

!  159. 

!  O.OOOE+00 

! 

O.OOOE+00 

! 

-1.17 

1 

! 

6 

-16.2 

!  -7.67 

!  150. 

!  O.OOOE+00 

i 

O.OOOE+OO 

1 

-5.35 

i 

6 

-0.416 

!  -8.09 

!  159. 

I  O.OOOE+00 

! 

O.OOOE+OO 

i 

-3.06 

1 

7 

4.04 

!  6.04 

!  43.0 

!  O.OOOE+00 

! 

O.OOOE+OO 

! 

-1.22 

! 

8 

-6.35 

!  -7.92 

!  169. 

I  O.OOOE+00 

! 

O.OOOE+OO 

! 

2.180E-16  I 

! 

0 

-7.49 

1  4.88 

!  38.2 

!  O.OOOE+00 

! 

O.OOOE+OO 

! 

-3.102E-16  1 

1 

10 

-7.49 

I  4.88 

I  38.2 

I  O.OOOE+00 

! 

O.OOOE+OO 

! 

-2.773E-16  ! 

i 

11 

-0.416 

!  -8.00 

1  169. 

1  O.OOOE+00 

! 

O.OOOE+OO 

I 

3.05 

! 

12 

4.04 

1  6.04 

1  43.0 

I  O.OOOE^OO 

I 

O.OOOE+OO 

I 

1.22 

! 

13 

21.8 

I  8.32 

I  60.7 

I  O.OOOE^OO 

! 

O.OOOE+OO 

I 

-1.03 

! 

14 

-6.74 

!  -6.77 

!  34.4 

!  O.OOOE-pOO 

1 

O.OOOE+OO 

1 

-6.16 

! 

15 

3.64 

1  3.34 

!  41.8 

1  O.OOOE-pOO 

! 

O.OOOE+OO 

1 

-7.60 

! 

16 

-3.09 

!  -1.62 

!  37.4 

I  O.OOOE-pOO 

1 

O.OOOE+OO 

1 

-0.46 
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J 

17 

!  -6.04 

!  -1.66 

! 

36.3 

! 

O.OOOE+00 

! 

O.OOOE+00 

! 

-6.74 

! 

! 

18 

!  -15.4 

1  -7.86 

! 

30.4 

! 

O.OOOE^OO 

1 

O.OOOE+00 

! 

-7.24 

! 

L 

19 

!  -10.9 

\  -6.67 

! 

32.9 

! 

O.OOOE+00 

! 

O.OOOE+00 

! 

-4.68 

1 

! 

20 

1  4.68 

t  4.18 

! 

42.6 

r 

O.OOOE+00 

I 

O.OOOE+00 

! 

-1.38 

I 

! 

21 

!  -21.8 

!  -7.27 

! 

159. 

! 

O.OOOE+00 

! 

O.OOOE+00 

! 

1.17 

1 

! 

22 

!  -19.8 

!  -8.08 

1 

159. 

! 

O.OOOE+00 

! 

O.OOOE+00 

! 

1.35 

! 

! 

23 

!  -26.4 

!  -3.32 

! 

160. 

! 

O.OOOE+00 

i 

O.OOOE+00 

! 

10.4 

! 

! 

24 

!  -16.2 

!  -7.67 

! 

159. 

! 

O.OOOE+00 

! 

O.OOOE+00 

! 

6.35 

! 

! 

25 

!  -10.9 

!  -6.67 

I 

32.9 

! 

O.OOOE+00 

! 

O.OOOE+00 

1 

4.68 

! 

j 

26 

!  4.68 

!  4.18 

! 

42.6 

! 

O.OOOE+00 

1 

O.OOOE+00 

! 

1.38 

! 

1 

27 

!  18.0 

!  2.20 

i 

46.9 

1 

O.OOOE+00 

f 

O.OOOE+00 

! 

-1.10 

1 

1 

28 

!  18.0 

!  2.20 

I 

46.9 

i 

O.OOOE+00 

! 

O.OOOE+00 

! 

1.10 

! 

! 

29 

!  -8.04 

!  -1.66 

J 

36.3 

1 

O.OOOE+00 

I 

O.OOOE+00 

1 

6.74 

! 

i 

30 

I-8.771E-02 

!  4.38 

! 

40.8 

1 

O.OOOE+00 

! 

0  OOOE+00 

-0.554 

! 

! 

31 

!  -3.45 

!  3.12 

! 

39.1 

! 

O.OOOE+00 

1 

O.OOOE+00 

! 

-6.081E-16 

! 

1 

32 

!  -7.21 

!  -9.31 

I 

159. 

! 

O.OOOE+00 

1 

O.OOOE+00 

! 

-1.40 

! 

j 

33 

!  -4.08 
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